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ABSTRACT 

Breakouts or preferred directional spalling along the walls 
of boreholes have been observed in many wells throughout the 
sedimentary basin in most parts of Alberta. Many studies 
have been carried out into the possible cause of this 
phenomenon. In this study, data focussing on the 
relationship of breakouts with depth, rock types and age of 
the rocks have been investigated in 50 wells scattered 
through most parts of Alberta plus one well in British 
Columbia. A regression of breakout azimuths on depths shows 
regression coefficients not significantly different from 
zero. Furthermore, it has been noted that breakouts have no 
relation with the rock types and the age of the rocks. The 
breakout azimuths in ae well are tightly grouped around a 
mean direction which, in most cases, 1S NW-SE. 

These observations and other evidences lead to the 
conclusion that breakout formations are better explained by 
the stress concentration at the walls of the boreholes, 
rather than the result of the oligos ial encountering 
pre-existing fractures or zones of weaknesses. From this 
Gonclusiony-= fhe «orventation: ©of “the smaller horizontal 
aoe nenyeall stress is NW-SE and the larger horizontal 
principal stress is NE-SW. The consistency of the breakout 
azimuths with depth, wherever these phenomena are observed, 
NayetMolyephatesuch large horizontal principal stresses are 
not necessarily limited to the top crustal sedimentary rocks 


but may be continued down into the deep and older 
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I. THE PHENOMENON OF ALIGNED BREAKOUTS AND ITS MEANING 
The primary purpose of borehole dipmeter surveys is to 
determine magnitudes and azimuths of dips of bedding planes 
from observations in a single hole. Prior to such _ surveys, 
the subsurface structural picture was determined primarily 
from seismic sections and correlation of marker horizons 
between three or more wells not drilled on a straight line. 
Unfortunately, geophysical surveys normally cover large 
areas and are good for determination of average regional 
dips butegive poor resolutvon of local dip variattons. In a 
like manner dip determination of marker horizons’ from 
correlation is valid only if the bed used is truly a plane 
Surface and continuous between wells. 

About 1970 the Four Arm Dipmeter (also known as_ The 
High Resolution Dipmeter, HDT) began to supersede the three 
arm type of tool. This thesis 1s largely concerned with 
Oil-wells of non-circular section. Figure 1 reveals that the 
HDT is capable of determining two different diameters of a 
non-circular hole and their azimuths, whereas the three arm 
instrument cannot do this. 

In the» *tour “arm dipmeter shown im Figure 2, the four 
arms are azimuthally spaced 90° apart. Each spring loaded 
arm carries a pad of electrodes which enables the electrical 
resistivity of the rock near the electrodes to be recorded. 
The four correlation curves in Figure 3 show typical 
resistivity variations which enable beds to _ be recognised 
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Figure 1.... Possible positions assumed by 3-arm caliper 
and 4-arm dual calipers in the same elongated hole. 
(schematic) 
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Figure 2.... Typical 4-arm dipmeter log record 
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values for a bed to give its dip in magnitude and direction. 

Resistivity changes on only one or two electrode pads may 

assist in identifying fractures as is illustrated 

schematically in figure 4. 

Traces not yet discussed in figure 3 are of principal 
interest in this study. On the far right the two orthogonal 
diameters indicated by opposed pairs of calipers are 
recorded. These provide information concerning the shape of 
the borehole. On the left we require the trace which shows 
the azimuth of No.1 electrode, relative to a magnetic 
compass housed in the tool. Other traces on the left are; 

1. The nearly vertical line which shows the deviation from 
verticality on a common scale from 0° to 9°, 

2. A dashed line which runs diagonally or rotates as _ the 
tool} 1s. drawn wp an stbhencaecular portions! of the hole. 
This line ceases to rotate on encountering an elliptical 
zone and gives the relative bearing of the tool ona 
scale of 0 to 360 degrees with respect to azimuth of 
No.1 elctrode. From this the bearing of the hole drift 
can be calculated. 

The tool normally rotates as it is drawn up the well and in 

many depth ranges the calipers will indicate essentially 

equal diameters equal to the drill-bit diameter. This can be 

Seen in the bottom and top sections ‘of - figure 3. It is 

observed, however, in other depth ranges that the tool 

ceases to rotate and the calipers indicate different 


diameters.  Thisewise ithe result “of “extensive (fracturing 
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cauSing the borehole to spall and wash out on opposite 
Sides,thus becoming elliptical. Usually the smaller diameter 
1s equal to the drill bit. This condition has been named by 
Babcock (1978) a BREAKOUT, and is illustrated in figure 5 D. 
Other situations encountered in a borehole are also shown in 
Bxguces . 5 VA to Cacinsfigure 5.0 one .pair -of ‘calipers has 
become trapped in a locally elongated section of the hole. 
The orientation of the breakout can be estimated from the 
azimuth of the trapped caliper pads. In Figure 3, a breakout 
is shown between depths 9418 and 9527 feet. Calipers 2 and 4 
(solid diameter line) indicate the greater diameter and the 
Stationary azimuth trace shows No.1 electrode near 200°. The 
AgomuUrnhe cols the yslong axis of. @his breakout @irs near 
200-90=110° relative “to the ‘magnetic. north. The final 
orientation is 110° plus the magnetic declination. 

The success of the HDT in determining sedimentary 
bedding planes has led to an extensive use of this tool in 
many areas, starting from the West Texas Ellenberger play in 
1969, the cretaceous of North LouiSiana and MissisSippi and 
also in the limestone of the Mooringsport of the Waveland 
field of Mississippi. 

During the interpretation of Schlumberger 4-arm 
dipmeter logs wells in Alberta, Cox (1970) observed that 
through the Wapiabi shale (a predominantly dark gray shale 
of tbetimbiGoliorado and post Colorado age ‘in the upper 
cretaceous period) interval of the Strachen-Ricinus area 


(see Fig. 16 shaded portion) of West Central Alberta 
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Figure 5.... Dual-caliper response to possible conditions 
of borehole shape. 
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adjacent to the Rocky Mountain foothills, where the borehole 
WasmenOtpwal obitwcizeJit=\ was ‘always elongated in. “the 
northwest-southeast direction. Structural dips were usually 
northeast or southwest. Nevertheless, several wells in this 
acecm ingewhich  (Gtructurals dips were difterently oriented 
showed the same northwest-southeast elongation. This 
Observation led Cox to make a broader statistical study of 
{arm aipmeter: Logs in, 317) zoness., from 17 «wells scattered 
about Alberta and the North West Territories. 

This new study ranged from Cretaceous’ shales EO 
Devonian carbonates, in which dip values ranged from 0° to 
30°. The orientation of the long axes of the holes had an 
average direction of N47°W to S47°E. A few elongations 
formed another group perpendicular to this major trend (Fig. 
Guay Thal the dip direction emsinot a controlling factor for 
the hole elongation is shown in figure 6 B. 

An interesting feature of these elongations observed by 
Cox (1970) was that they were confined to discrete depth 
intervals within a hole and separated by uncaved beds. This 
observation led Babcock (1978) to suggest that the 
elongation might be caused by the drill encountering zones 
of steeply dipping fractures. Babcock later undertook a 
Study of the phenomenon of hole elongation using the caliper 
portion of the 4-arm dipmeter logs, and extending his 
investigation to 23 wells in various sedimentary rocks 


widely distributed through Alberta. 
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Figure 6.... Orientation of breakout long axis with respect 


to (CA) North and *(B) Dip azimuth. (J. Cox, 


1970) 
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He noted the following properties of breakouts 

1. The breakout zones may be short discrete events, or may 
persist over depth intervals of several tens or hundreds 
of feet. 

2. A breakout of 2 to 4 feet commonly slows down the tool 
rotation or may have no noticeable effect on rotation; 
whereas a longer breakout 1S commonly associated with a 
cessation in the rotation of the tool. 

To explain his observations, which were consistent with 

€ox'sa(1970) observations,’ Babcock! turned to» his’ studyo of 

the regional jointing in exposed bedrock of the Alberta 
plains and in the Fort McMurray area of northeast Alberta 

(Babcock, 1973,1974,1975). Within this area two orthogonal 

joint systems made up of vertically dipping extensional 

fractures are present. System 1 has sets striking northeast 
and northwest roughly normal to and parallel with the Rocky 

Mountain Belt. The northeast striking set iS dominant in the 

sense that it 1S most commonly developed in that joints of 

this set usually cut across joints of other sets, which 

terminate against them. A second joint system (System 2, 

Babcock, 1973) having sets striking roughly North-South and 

East-West was also observed. This system was common in_ the 

Fort McMurray area and in the southernmost part of Alberta 

and occurred elsewhere in the province. These two joint 

systems were reported (Babcock, 1973) in nearly flat-lying 
sedimentary rocks ranging from late Devonian to Paleocene in 


age, in outcrops of shales, siltstone, sandstone, limestone 
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dolomite and coal distributed through many thousands of 
square kilometers and several thousands of meters of the 
Stet doen ec, eCOlLUMn A OUtm Of al vtotal of 4 1)236) nearly 
vertically dipping joints measured, Babcock observed that 
although joint sets were well defined at most stations, set 
directions between stations varied considerably. 
Nevertheless, the sets striking northeast and northwest were 
most developed. 

Regions covered by Babcock in his study of breakouts 
included the Claresholm district. Here he compared the joint 
azimuths measured at 10 outcrops and the azimuths of 
elongations from five wells. The mean directions of the 
elongations he observed showed a variation in strike ranging 
trom, 129° (‘to -@5505",0 thuseecorresponding: closely’ to the 
northwest-southeast striking joint set. However, the 
elongations showed much less variability than do the strikes 
of individual joints of the northwest striking set in the 
area. A comparison of all azimuths of elongation (41 in all) 
with that of all joint strikes (1,236) in the area is shown 
in Figure 7. Clearly there is no close correspondence as a 
result of the predominance of one joint set striking at 
approximately 175° (North-South), which falls under 
Babcock's System 2. 

Similar comparisons were made by Babcock in his Red 
Deer study between 400 azimuths of joints in bedrocks and 67 
breakout elongations in nearby oil wells. Here, however, 


joint sets belong to the System 1 and the dominant direction 
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Figure 7.... Histogram comparing (left) strikes of 1,236 
joints measured at 10 outcrops in Claresholm region, 
southwest Alberta, with (right) azimuths of 41 elongated 
breakouts from five wells. (Babcock, 1978) 
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of hole elongation was northwest-southeast at about 140 °. 
Babcock's study of 23 wells scattered throughout 

Alberta led him to draw the following conclusions,some of 

which had been made by Cox (1970): 

1.” “Ine most sof Grthe 23. wells istudiedpcthe ‘breakouts ¢ane fin 
the northwesterly direction. However, in some wells a 
small number of breakouts may follow a northeast trend. 

2.) etthe azimuths of elongation are parallel to well 
developed sets of regional joints which are present 
wherever joint directions have been studied. 

3. Hole deviations from verticality, changes in the 
magnitudes and directions of bedding dips, and lithology 
are unrelated to breakouts. 

4, The preferred orientation of breakouts could not be the 
result of the drill encountering vugs because this would 
lead to hole elongations with no preferred orientation. 

5. Babcock further Stated that through his informal 
discussion with engineers involved with hydraulic 
fracturing and from consideration of the tectonic 
setting, the inferred direction of intermediate stress 
in rocks of the Alberta plains is normal to the trend of 
the Rocky Mountain Belt or roughly Northeast Southwest 
and the maximum stress is vertical, with the minimum 
stress direction inferred to be Northwest Southeast 
parallel with the dominant azimuth of hole elongation. 

From these general conclusions Babcock finally concluded 
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encountering steeply dipping fractures or zones of steeply 
dipping fractures, which may or may not be open. The 
oversize holes in carbonate rocks he related to 
solution-widened joints. 

Bell and Gough (1979) found Babcock's hypothesis for 
the cause of breakouts very difficult to accept. Their most 
serious objection was the fact that the azimuthal 
distribution of breakouts is very unlike that of joint sets 
at the surface. They argued that given four concentrations 
of joint directions on the surface (Babcock 1973) NW, NE, N, 
and E, they would expect concentrations of breakout azimuths 
in all the four directions on Babcock's hypothesis, rather 
than only one significant concentration of subsurface 
breakout direction (NW). 

Secondly, Bell and Gough wondered if such surface joint 
directions would necessarily be parallel to those at depths 
extending beyond 2 km at which breakouts have been observed. 
In this connection they remarked that sediments involved in 
the study in the Alberta Plains were laid down in varying 
tectonic settings between Devonian and Cretaceous times; and 
furthermore joints may arise in several ways after 
deposition of sedimentary rocks. 

In an attempt to find an interpretation that would 
account for both the breakouts in the wells and the’ surface 
joints System 1 of Babcock, Bell and Gough (1979) proposed 
that both features result from a general stress field acting 


throughout the Alberta Plains and oriented with the larger 
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principal horizontal stress NE-SW and the smaller of the two 
horizontal stresses parallel to the Rocky Mountain fold 
axes. Using the Kirsch equations for the stress near a 
Circular hole, in a medium under biaxial stress, they argued 
that with large enough, unequal horizontal stresses oriented 
as above, the holes themselves could concentrate the 
stresses Sor asi’ tol’ produce®subsurface breakouts¢with*long 
axes in the NW-SE direction. Such concentrations of stress 
are believed to arise at the time of drilling. 

If the breakouts in Alberta are shear fractures caused 
by the concentration of stress at the walls of the borehole 
as proposed by Bell and Gough (1979), the larger horizontal 
compression is orthogonal to the azimuth of the breakouts 
and hence nearly Northeast-Southwest in Alberta. Their 
hypothesis provides no indication of the magnitudes of the 
principal stresses other than the fact that the horizontal 
principal stresses are large and unequal. Gough and Bell 
(1982) have recently made a quantitative study of shear 
fracturing near a borehole, and have shown that such 
fractures can form breakouts with the observed properties. 

Stress measurements in mines in several continents 
reveal that the vertical principal stress approximates 
(within 20% or less in most cases) the overburden pressure 
(McGarr and Gay, 1978). Furthermore, the proximity of the 
thrust faults of the Rockies to areas covered in the studies 
by Cox (1970) and Babcock (1978) coupled with the fact that 


the Rockies are fold mountains led Bell and Gough (1979) to 
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suggest that a situation with S, (the maximum principal 
SGpess) iverticality avenormab+faults stress titeld)wlis: “most 
unlikely in Alberta. 

This fact is further strengthened by the occurence of 
breakouts at depths of 500 meters or less. Laboratory 
experiments by Gay (1973) indicate that the threshold 
pressure for spalling in a circular hole under equal 
orthogonal pressures is about 100 MPa. With unequal 
pressures, however, spalling would be expected with the 
larger compression around 30 MPa. But at a depth of 500 mm 
the overburden pressure will be about 12 MPa. If S, is 
vertical then both S, and §,, which are horizontal, will be 
smaller than this value (12 MPa) and may be too small to 
produce shear fractures. This now leaves uS with either a 
thrust fault or a strike-slip (wrench) fault stress field. 
Bell and Gough's hypothesis fits either of these two. stress 
fields and the breakouts are consistent with either case.It 
is worth noting that the breakouts alone cannot reveal which 
of the two stress situations prevails in Alberta. 

In a later paper, Gough and Bell (1981) added 
additional wells in Alberta to their study of breakouts.They 
also extended the area covered into northeastern British 
Columbiayovinroadditionsresults«fromehydraulic fracturing in 
an oil field in West-Central Alberta were discussed. They 
suggested that some local anomalies in breakout orientations 


might represent local variations in the stress field. 
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Eonand -Morton,..(1976)» attributed. failures. -of.. ‘tunnel 
BOOoLSs 1n)eOntanio to interaction of the tunnel with large 
horizontal compressive stresses transverse to the tunnel. 
Dusseault (1977) suggested that unequal horizontal 
compressive stresses prevail in several areas in Alberta, to 
account for the propagation of vertical fractures during 
hydraulic fracturing in Northeast-Southwest vertical planes. 
Such fractures imply a horizontal least principal stress 
(S;)directed Northwest-Southeast. The formation of tensile 
fractures normal to S, during hydraulic fracturing was first 
proposed by Hubbert and Willis (1957). This leads to the 
expectation that in .a thrust -stress.field (S, vertical) 
hydraulically formed fractures will be horizontal (Kehle, 
1964). This is difficult to verify in many cases, because 
horizontal fractures are hard to detect. Zoback et al (1977) 
argue that the inflatable packers inhibit such horizontal 
Bractunes, sandethatehor zontal fractures, formseonly aif» fluid 
penetrates pre-existing planes of weakness. This contention 
is supported by laboratory experiments by Haimson and 
Fairhurst (1970). Haimson (1976b) suggested that if S,= S; 
then the induced fractures initiate in a vertical plane and 
then become horizontal as they propagate based on the fact 
that energy propagates along the path of least resistance. 
This can be observed in the pressure-time history where the 
asymptotic Shutsin pressure is smaller than the 
Lnotankaneous sShuts ine preSSure mw. line-suchys.a Situation ithe 


instantaneous shut-in pressure is not equal to the least 
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principal stress but rather the smaller horizontal principal 
Stress with the asymptotic shut-in pressure being equal to 
the least principal stress S,. 

In the J Lease of the Pembina Oilfield, Macleod (1977) 
noted during secondary oil recovery from Cardium Formation 
sandstone that permeability was greatest in the 
Northeast-Southwest direction. This observation implied 
either a fracture system or else a permeability trend 
oriented Northeast-Southwest. But Neilson (1957) showed the 
isopach axes in this formation and area to be trending 
Northwest-Southeast. Thus one would expect a permeability 
trend related to the process of deposition to run NW-SE. 
Gough and Bell (1981) therefore concluded that Macleod's 
observation could be better explained in terms of vertical 
fractures oriented NE-SW, which would imply a NW-SE 
Orientation ofpeS,a!'<Six wells loggedswith the: HDT »inithis 
area were identified and studied by Gough and Bell. The 
breakoutsS occur in rocks of Paleozoic ages over depths from 
2500 to 4500 meters. Azimuths of hole elongation were 
tightly grouped with mean values ranging between 132° and 
and 141°, with variations over approximately 208 in 
individual wells. 

Imperial Oil in 1978 published information on _ their 
Cold Lake projects in far eastern Alberta. Induced vertical 
fracture orrentations varied from N 30° E to N 45° oF 
implying that the minimum principal stress was horizontal in 


this region and oriented between 120° and 135°. The larger 
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horizontal stress would then lie between 30° and 45°. 
iinipemal Oil “ran no HDT logs in wells in the’ Cold Lake 
region (J.S.Bell, personal communication with the Company). 
However, Gough and Bell (1981) located one well in this 
region logged by the HDT. The azimuth of a single breakout, 
at depth 477 to 479 m in the Clearwater Formation, was 131° 
Suggesting the major horizontal stress direction of 41°. 
Recently, Gough and Bell (1982) have reported comparisons of 
breakouter azimuthsaumns Coloradost innt i Basta Texase@and in 
northwestern Canada with other indications of stress 
Orientations in those regions. In every case the 
observations support their explanation of breakouts as a 
consequence of stress concentration by the borehole in an 
anisotropic stress field. The hydraulically induced 
fractures at West Pembina and at Cold Lake indicate a 
Strake-sino gicypee istiressoe fieldsea(Ss esandkSsehorizontady Sh 
vertical) is most likely present in most parts of Alberta. 
This thesis reports a study of the statistics of the 
azimuthal distribution of breakouts in oil-wells in Alberta 
in relatmon tordepth, Wuthology; formationrandrage ofsrocks; 


together with the geophysical implications of breakouts. 
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II. STRESS AND FRACTURE IN THE EARTH'S CRUST 
STRESS 
If one considers a simple prism of cross-sectional area 
a subjected to a force F,as shown in Fig. 8, the stress S, 


acting on the end surface ABCD is given by 


It can be seen that the force F hasS no component acting 
parallel to the ABCD surface. This means it exerts no- shear 
stress on this surface. By definition, any stress acting 
perpendicular to a surface along which the shear stress is 
Zero 15 a’ principal stress. Thus. in this case Sz is a 
principal stress. Here, the suffix indicates the direction 
in which the stress acts. Other principal stresses may be 
orientated parallel to the X- and Y- axes, and would be 
designated Sy and Syowhespectively. If the relative 
intensities of the principal stresses are known, they may be 
termed the maximum (or greatest), intermediate and minimum 
(or least) principal stresses (ives Soa, Sey Sz 
respectively). 
If one considers the action of the force F on a surface 

GHIJ inclined at an angle @ as shown in the Figure 8, the 
component of F acting normal to GHIJ is given by, 

F, = F.sing. 
However, it will be seen that the area a'of GHIJ is. greater 
than the area a of ABCD, and that 

a’ = a/singd<. 


Therefore, the normal stress S, acting on the inclined 
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Figure 8.... Normal and Shear stresses acting on external 
and internal surfaces of a unit cube subjected to a force F 


a3 


Surface is 

FEF, /a’ = F.sin@/a' =(F/a).sin?@ 

so that 

Ss “EWS, Csinte jt-etese reste eres ——- Rn -- eb ee 2-2 

Similarly the component of F tangential to the inclined 

plane is given by 

F, = F.cos@. 

Consequently the shear stress (T) acting along the plane 

equals 

F /a' =F.cos@/a'= (F/a).cos@.sing 

that is 

RE COS SL 2-3 
Shear stresses as well as other stresses may be related 

to co-ordinate axes. We denote by T,, the shear stress on 

¥eqolanes), acting ~ parallel .'te *the "y-axis. “Similarly T,, 

denotes a shear stress acting in the zx-planes parallel to 

the x-axis. The stress across a plane whose normal is in the 

yecirec tions sa Wiihave. ‘components T'S) = Ty, and that 

across a plane whose normal is in the x-direction will have 


components rs jo Tyo Ty LOThe nine quantities” below 


xy XZ 
Sy Txy Tz 
Tyx Sy Tyz ce = a ee Oe ee ee age ea ae ee 2-4 
og) ee & S; 


are called the stress components and give a complete 
specification of the stress at a point. It can be shown that 
Desay > 6Ty:6=Ty rand Ty HT"; SO that “only “six -quantities 


out of the nine are needed to specify the stress at a point. 
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The stress Component | Olu phoma +4), are avin facts. the 
components of a mathematical entity called a tensor,(and 
tensor analysis is much used in developing the higher parts 
of ‘the ‘theory of elasticity where’) x,y,z, 9 are. xreplaced.. by 
1,2,3, respectively). A tensor in which Beis 
TT. 9=T,, “ls said tobe esymmetric. 

In a biaxial stress field it can be inferred from Eqn. 
(2-2) and Eqn. (2-3) that the stresses acting normally and 
tangentially to a plane inclined at @ to S, and hence (90+@) 


to S,; will be represented (on superposition) by 


SS, 2sin?g+S3se0s*? 
and 
Toei (SLR Sane cOse~ ~ ee enee ss -Tess ooe ERS Se caeeel 2-5 


which can be written in- terms of the double angle 2¢ as 

S, ae SatSu0y 2h (NSeeS Jeicos2f)/72-and 

Peete Sono Oey 2a EEO Ee ea ee 2=6 
Equation 2-6 may be represented graphically by means of 
Mohr's stress circle, which will be discussed later in this 
chapter. 

Subsurface rocks are normally in a state of compressive 
stress except very near the surface because of the weight of 
overlying rocks. This overburden weight creates stresses in 
both the vertical and horizontal directions. The stress 
field at a point can be represented by S,, S2, S; as already 
defined. Over long periods of geologic time the earth has 
exhibited an appreciable degree of mobility during which 


rocks have been repeatedly stressed to:the limit of failure 
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to produce faulting and folding. The problem of rock squeeze 
(i.e. continuous deformation of the rock) pee been 
Gecogntsed bain gpantstmhofetthes cuustyajIninSotthern Ontario 
structures built in rocks have exhibited distress in various 
degrees during or subsequent to construction (Lo, 1978). 
Geological processes such as faulting,folding, and pop-ups 
or buckling of the surface rock Strata, without any apparent 
change of external loading, may be interpreted as evidence 
of existence of high horizontal stresses. In a study of 
geological features and movements of structures in rock in 
New York State, Rose (1951) associated the rock squeeze with 
the possible existence of high horizontal in-situ stress in 
the rock formations. Such phenomena require that substantial 
differences must exist between the principal stresses. Sbar 
and Sykes (1973) have used information from postglacial 
geologic features such as eaSt-west rock squeeze in western 
New York and Niagara, and bridge abutments moving together 
to infer that the maximum compressive stress in this area 
has an easterly trend. Postglacial buckles or pop-ups near 
Chippewa Bay, New York have also been ascribed by Sbar and 
Sykes aS not due to any environmental factors other than 
large horizontal compressive stresses. Pop-ups ~have also 
been encountered elsewhere in western New York where the 
lithostatic load has been reduced by quarrying (Sbar and 
Sykes personal commun. with J. Davies, 1972). Deformation 
indicative of a high horizontal compressive stress was also 


observed in various parts of Ontario (Coates, 1964), where a 
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pop-up about ctwoh andwhone thalfirmeters!{highye striking 
northwest, occurred overnight in a quarry after it had been 
excavated: to cac\ depth cofy, 155 meters. Finally the results 
obtained from geologic features, in-situ stress measurements 
and earthquake distributions in eastern and central North 
America may be related to the presence of high stress. 

The orientation of the trajectories of the principal 
stresses in space is largely determined by the condition 
which they must satisfy at the earth's surface. This is a 
free surface, on which normal and shear stresses vanish. 
Since the only planes on which the shear stresses are zero 
are those perpendicular to the principal stresses, it 
follows that one of the three trajectories of principal 
stress must end perpendicular to the surface of the ground, 
and the other two must be parallel to this’ surface. 
Therefore, in regions of gentle topography with simple 
geologic structures the principal stresses should be nearly 
horizontal and vertical, with the vertical stress 
approximately equal to the pressure of overlying material. 
Measurements of stress in mines in several continents’ show 
that in many cases the vertical normal stress is, within 
about! 20%; .equaletostheioverburden; pressure: Syo =): psg.h ., 
where p is the density of the material, g is gravity and h 
issthe depth. If the lateral, stress were due solely. to 
elastic: i response) of», the* rock-Sy,s thei horizontad, principal 
stresisesrywould \bese S=s= ov/(1-1)2.Sye =: Sy/ 3g wheres. or (AS 


Poisson's ratio. Such a stress field would produce normal 
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faults at any depth at which the stress difference exceeded 
whe -strengthk LOrmathelerock® However, observations in 
underground structures, quarries and hydraulic fracturing in 
boreholes indicate that the horizontal stresses S and s S>s 
are not generally equal and do not vanish close to the 
surface, a necessary consequence if the horizontal stresses 
were duesmsimplyi eto Pelastweysresponse ofS theeerock mito 
overburden pressure. 

MeaSurements in several continents of in-Situ stress in 
rocks show that in many cases the horizontal stresses are 
much higher than the vertical stress calculated from the 
overlying weight of the rock (McGarr and Gay, 1978 ; Hast, 
1973; Herget, 1974). Frequently the horizontal stresses § 
and s increase linearly with depth. It is thus evident that 
large and unequal horizontal principal stresses are 
widespread in the continental crust of the earth. 
Relationships Between Principal Stresses at Failure 

Most substances behave elastically at low stresses. As 
the stress is increased the body will begin to yield at some 
Polnteliarc, m@srauctwle Swhitemitethe@bodyeis brittletatiwidd 
fracture at some point without appreciable yielding. The 
term fracture implies the appearance of distinct surfaces of 
Separation fin Sthe!, body andoyield® is®usedtor) the>onset: of 
plastic deformation. Flow is used for unrestricted plastic 
deformation. Both flow and fracture (also known as failure) 
are observed on large and small scales in geologic material 
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interpretation of geologic phenomena. 

For brittle materials in tension, tensile or cleavage 
fracture takes place across a surface perpendicular to the 
Givecuron Of LenslOonai bigs JOA). sFor brittle “materials: ‘iin 
compression, shear fracture takes place along a pair of 
planes or a cone approximately in the direction of the 
greatest shear stress but always between this direction and 
the direction of the largest compressive stress ( Fig. 9 B). 
The maximum shear stress theory, which dates back to Coulomb 
Vi 3) sretessthat=fallure occurs at a \point when? the 
maximum shear stress is equal to the shear strength of the 
material.) 1f.S,2S,2S,, are the principal stresses at a point, 
taking compressive stress positive, the maximum shear stress 
has) magnitude §=(S,-S,)/2 , and occurs across a plane 
containing the direction of S, and whose normal bisects the 
angle between the greatest and least principal stresses. The 
theory implies that if C, is the compressive strength of the 
material in an unconfined compression test in which S;=S,=0, 
Se=Co pe thesematetial will fail across any’ plane inclined at 
To -SrOutne direction ~ol- compression- jfor- a) frictionless 
material. For real materials the angle is less than 45°. The 
theory also implies that the tensile and compressive 
strengths are equal. The theory has been modified by Navier 
to fit qualitatively most of the observed facts, and is used 
by “Anderson in a discussion of the types of geological 
fAuLeing ew liars moditication, Winstead «of assuming wethac 


fracture takes place across the plane over which the shear 
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Figure 9.... Failure of a brittle material.(A) Tensile 
fracture and (B) Shear fracture 
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Stress first becomes equal to the shear strength S, of the 
medium, it assumes that this shear strength is increased by 
a constant u times the normal pressure across the _ plane. 
Because of the analogy to ordinary friction in which the 
tangential force is u times the normal reaction, u is known 
asmprnesicoecliscnentectlinternabnirretion. Thus, Tf SplandhT 
are the normal and shear stresses across a plane, fracture 
takes place on the plane at which the magnitude of T first 


becomesTequaletoyS prtducss. thherefore; 


This 1s the Navier-Coulomb criterion for shear failure. Let 
us discuss the case of two dimensions first. If S,, and S, 
are the principal stresses,we can then write Eqn. (2-6) as 

SA =PhS. +52) 72SttokS, oS. kcos2gA2 

Re=pS noes esenee 7/2 tre fer See Sse Se --sS- 220 
Only values of @ between 0 and 90° need be considered, since 
eoplyurheimagmbude of bri toccursackmreEqmiti(2=7), psenmethat 
changing the sign of @ only changes the sign of sin2¢ and 
deestnot at bect tthe smagnitude sof »S,\ oomph. iThetishear Ustréss 
at failure is thus symmetrical about @=0. It can be shown 
that failure takes place across a plane whose normal makes @ 
Wath reSchewhereevtan2d | 4=i/u. tht weOry thene@=457eariocust; 
20=135 Yeand<O=67.52> thE Baca iinfiimity ;@o-> o0mthat rscas neu 
increases the plane of fracture moves towards the direction 
of maximumestressscValbues of cureofi athe ‘corder. of lis are 
inferred from the directions of fracture of rocks in testing 


machines, and also from geological faulting (Jaeger,1956; 


oy 

vd boeastonil at Asetinaae sheds 
ye 

,eneiq eds 266128, swensiq” ‘ti 
aij foidw ai noobs s5iva. grasibie Poa pee S 
: 2 Soe 

nwonk et o ~noiyose7 Lammiion ota aork 


T Bae ,2 22 , eu 2 aondoe fanssant 40-4 
a) iy 7 = 
vgos 33 ,ensig 6 ao I]8 eeeanite reede tan. 


> mA i 
2712 T 30 sbusingem edz Goidw Je ‘onal od 
3 AL : 41/3 Ak a 1; 
; * a ee 
a0 igisienT . gavel * 

: we eise 7 

f : 7 - 22) : a i 
| an ea we were et ee eee ee oO _ ee ee a uae 
ja? .esyuitea? seeda 163 noisesine énotueo~sebvalt sit 


serth snoienenié owe > deeo of7 - “eet 


: i : ae 
» ie Pb he W iSa42 = iad 4, poage: 332 sms af 
/ ee +. - ei 
S\Gsacok . 2.2) * NM git 


* 


a gan nai tk sya -¢\nsiie! aes 
inca *ye bail. asewsed, by 6 a re 


Of. (¥ -¢} 
sxga: segca, oft .T mn ye io usingen. ere 
> ao f . (7 a hs 4 5 
agdem Lamon. seortu maannatl —oe sa 
pra wees 


noizoe1ib srg abyewo3 29von ents 7 
Ms 


“a 


cs DE, Ohee nedd Onu aI- 


oe @e@ et pais ge<--8 suuiazank 


oe ft vebse: ott toe 


Sal 


Price,1973). It should be remembered that,because of the 
symmetry in @ mentioned earlier, the theory leads to two 
possible planes of fracture equally inclined to the 
principal stresses and gives no reason for. preferring 
either. For failure under combined stress in terms of the 
compressive and tenstde strengths, CC. and T, of the 
material,the relation 
Cy P= Ne ea) (fcaeet)) =u) 
holds. The theory predicts that the compressive strength of 
a material is always greater than its tensile strength but 
the ratio 1s rather smaller than that found in practice. The 
theory also predicts that under any conditions the normal to 
the plane of fracture makes the same angle arctan(1/u)/2, 
with the direction of greatest principal stress. Although 
this 1S approximately true for compressive stresses, it 15 
very far from the truth in the case of pure tension when the 
failure is uSually by brittle fracture with the plane of 
fracture normal to the direction of tension. The reason may 
be that T, should not be the actual tensile strength but the 
value at which shear failure in tension would take place if 
tensile fracture did not occur in practice before this value 
is reached. Nevertheless, the theory gives a reasonably 
accurate account of the behaviour of rocks under combined 
compresSive stresses. 

Mohr's representation of stress and failure is useful 
(Fig. 10). Mohr's theory assumes that at failure across a 


plane the normal and shear stresses across the plane, S, and 
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Figure 10.... Representation of stresses on any plane ina 
two-dimensional stress system by means of Mohr's stress 
circle. 
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T, are connected by some relation T = £(S,). This. relation 
may be plotted on a (S,, T) plane and since changing the 
Sign of T simply changes the direction of failure but not 
the condition for it,the curve is symmetrical about the 
S,7~axis. Any state of stress can be represented by a Mohr 
circle on the (S,, T) plane. In representing the normal and 
shear stress across a plane in a biaxial stress field on the 
Mohr's Stress circle, the normal stress S, and shear stress 
(T) are chosen as co-ordinate axes. Here the greatest 
principal stress S, and the least principal stress S; are 
represen ede by min and™OM respectively on *'the- oS? -axis. "The 
quantity (S,+S,)/2 of Eqn. (2-6) represents the mid-point C 
on the S, axis between M and N, while (S,-S;)/2 represents 
half the distance between M and N. If a circle is drawn with 
centre C and radius (S,-S,;)/2 , then for any specific values 
ers, ®and SS.) *5thi's circle’ “represents ‘the ‘condityons’ of 
Equation 2-6 where 2% is measured as indicated in Fig. 10 a. 
This construction as we shall see is used in representing 
the values of shear and normal stresses at failure. If this 
circle lies wholly within the failure envelope ABA'B' in 
Fig. 10 b, the stresses involved nowhere attain the critical 
values. If any portion lies outside it the material could 
not withstand the stresses. Obviously, the limiting case is 
that’ Of a’circle such as those of centres C,, O and C, which 
just touch the curves AB and A'B'. In this case failure will 
take place under stress conditions corresponding to the 


points PP', that is, over planes whose normals are _ inclined 
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at angles of half the angle PCN (Fig. 10a) to the direction 
of the greatest principal stress. The clirve AB will be the 
envelope of all the circles corresponding to all conditions 
at which fracture takes place and for this reason is known 
as the Mohr envelope. 

In principle, three circles which touch the envelope 
can be found from simple experiments. These are those of 
Centres C;, O’and C, (Fig. 10 b)”, corresponding to tension, 
Simple shear, and compression. But in practice it is 
Gr fhpceult to perform shear or tensile tests on_ rock 
materials. An approximation to the Mohr envelope for many 
rocks, *"rrom-> results °ot -*triakial “ tests, °-rs' the pair o¢€ 
straight lines 

Pg) Seep 
corresponding to the Coulomb-Navier criterion for fracture 
under shear stress, where T, is the shear strength. In this 
case the normal to the plane of fracture makes an angle (in 
the second quadrant) arctan(1/u)/2 with the direction of the 
greatest principal stress. When applied to triaxial stress, 
the Mohr theory leads to the result that only the Mohr 
Circles’ for the plane “containing ‘the greatest “and least 
principal stresses need be considered and that fracture 
always takes place in planes containing the direction of the 
intermediate princroal’'"stress’.. *9This »“is'not* ‘eleogether 
consistent with the experimental results. 

One important application of these results is the study 


of faults which are fractures of the rocks of the earth's 
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crust (Fig. 11). Geologists distinguish three major types of 
faults and Anderson (1951) has shown that these are 
determined by the relative magnitudes of the principal 
stresses. One principal stress will always be vertical, and 
three caseS arise according as this is the greatest, 
intermediate, or least principal stress. In the case where 
the vertical principal stress is the least in magnitude and 
the other two principal stresses are compressive (Fig. 11 
A), the fault is known as a thrust fault. The planes of 
fracture pass through the direction of the intermediate 
principal stress and makes angles of less than 45° with the 
direction of the greatest compressive stress which is 
hoprzconva Wee (hig. |) Bes Thissapplies in regions “of “active 
tectonic compression or erosion or vertical unloading. When 
the intermediate principal stress 1S vertical, we have 
Figure 11 C. The failure can take place on either of two 
vertical planes (Fig. 11D), which are equally inclined at 
anglesmeot es less.  thanry 45° towthe direction on the i qreatese 
compressive stress. This type of fault is called a 
transcurrent, wrench or strike-slip fault. When the vertical 
principal stress is greatest as may often be the case at 
considerable deptns,) ands ity srimemctructunressenear ge the 
surface, the type of failure is called normal faulting. Here 
the faults make angles less than 45° with the vertical (Fig. 
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Figure 11.... Types of faulting in relation to the 
principal stresses; Avand B, Thrust fault, -C and D;-Wrench 
fault, E and F, Normal fault. 
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Stress Concentration Caused by the Borehole 

the-s presence of] a borehole distorts the pre-existing 
stress field, because as discontinuities in the earth's 
crust, boreholes cause redistribution of stresses in the 
nearby rock. An approximate calculation of this distortion 
is made by assuming that the rock is elastic, the borehole 
smooth and cylindrical and the borehole axis vertical and 
parallel to one of the regional principal stresses. Kirsch 
(1898) found an analytic solution of the problem of the 
stress field near a small hole of radius a ina large plate 
under uniaxial compression S. The solution as given by 
Timoshenko and Goodier (1951) is given below. 

Let Figure 12 represent a plate subjected to a uniform 
compression of magnitude S in the x-direction. If a small 
circular hole is made in the middle of the plate, the stress 
Seen Nee in the immediate vicinity of the hole will be 
changed. However, by Saint-Venant's principle we know that 
the change will be negligible at distances which are large 
compared with a, the radius of the hole. Let us consider the 
portion of the plate within a concentric circle of radius b, 
large in comparison with a. The stresses at the radius b are 
effectively the same as in the plate without the hole and 
are therefore given from Equations (2-2) and (2-3) by 
( Sp eb = S-cosS?M= S/2.(1+c0s29) 

WM ee eo ANIC ate aa eer een ag 2-9 
For Stresses disposed symmetrically about a point of 


weakness in the earth's crust the direction of the principal 
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Figure 12.... A plate submitted to a uniform compression of 
Magnitude S in the X-direction with a small hole in it 
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Stresses will be vertical, radial, and tangential. The 
forces acting around the outside of the ring in Equation 
2-9, having the inner and outer radii r=a and r=b give a 
Stress distribution within the ring which can be considered 
as consisting of two parts. (a) The first is due to the 
constant 1/2 S of the normal forces. (b) The remaining part 
consists of the normal forces 1/2 S.cos2%g together with the 
shearing forces -1/2 S.sin2%. The stresses produced by this 


second part may be derived from a stress function of the 


form 
Tee OME ce Penne hae Oe SO ee ee ee 2-10 
This stress  functyron can be substituted into the 


compatibility equation 
(+1241 2,)( Bes 120 + 4FF)-0 

dr2 or Orr? SGA ror re 3d 
to yield an ordinary differential equation from which f(r) 
can be determined. 

(Se? Jor *td/rar-4/r*) (a*t/dar?+ at/rdr-4f/r? }=0 
The general solution is 

Mayes Agel FF Bore + Cee FD 

Equation 2-10 then becomes 
etme eee eC, i FD) cOSe eee ee al Oe 


The stress components are; 


2 
Settee se POD C 4D 2¢ 
ama (2A + or + 45 ) cosine 
2 
S SS om = 2 6C (ois eee aes cles "ae es asa oe ps aa ool eo 2-12 
Eis ( 2A + 12Br + C, cosine 2¢ 
oe ee + 6Br2-6C - 2D )\sine2¢ 
? Ara 2s ags eae) 


ee 

ad? + fetan neonee ‘tn waniiiee: 
; _ hi staal 

& ave ds; bre oat bins pelea 

bersbienos ed nas tobi ie pata ‘of nike 


otyeupa wk aatt. 


bei 


hp 


ei3 07 aub ef: 289 £2 oh? (a) : efveq a, 


‘adh 
2 


Shintames sdf (ad) 2 99703 Learn oe 99: 
a phe 


gta daly lonsepos Oires.8 SY | sore’ keirxon 


ar ‘ Db 3) 3 a 
aid yd Gsoubody egaretee ad fina: e\t+ 
; . eMC (eles = mA Le 
aia io notsseuw? as e172 2° Ox bevizeb Cdl yom” 
fa Li : - 
a : (mg 
oe = ie: ay o 


Py : ; me 6 
= => es ee ee ee ee 1 toon ta) bet 


+e 
; ao | pore a ‘aaouse on 


i? 
ie 
Seb 
nd 
he 
, 
a 
fu 
uw 


es Pragins sibiat am 


~ - a - 
( ; ‘ez: Pos Oe Se Sed “+ & eo i s) 
ds & sy ve i % . S&S § ' . oe he ee -— re 
: aii Z) . i 7 * ; fia : | “ Wy 
ityt dotdw. movt mor doeutie isetnars2 216 exsetbrs all bie ied Xs 


bonimrsseb 


Ge (* ¥\eb-obs\ eb eae i a oe] f bi pyar abe 38 ff 


ce Ti het nutos isxene 
\ ch 7 an ; 


A * ce + a oe i: 


ee ee 


40 


The constants of integration A,B,C,D are determined from the 
conditions of Eqn. 2-9 for the outer boundary and from _ the 
condition that the wall of the hole is free from external 


forces. These conditions give 


Pret OC Ds + AD/ be = -S/ 2 

2A + 6C/a‘* + 4D/a? = 0 

eho OBO? —"6C/b* = 2D/b* = S/2 
2h + 68a°. — BC/a* — 2D/a? = 0 


For an infinitely large plate a/b = 0. These four equations 
can then be solved to yield A=-S/4, B=0, Ca— aes (A: 
D=a?.S/2. These constants can now be substituted into 
equation 2-12, and by adding the stresses produced by the 
uniform compression S/2 on the outer boundary we obtain 

SES oe ee 5/2 Sate * = Aan ).cosed 
aS radial stress 

BS ea ata aS aa COS 2 a ge 2-13 
as the tangential stress and 

eo Sen lo3o tir (49a er 7a oS 1 nee 
as the shearing stress, where the stresses are expressed in 
polar co-ordinates with the centre of the hole as the origin 
and the plane stress components are at a point (@, r) 
exterior to the hole of radius a in a plate under uniform 
uniaxial stress, with @ measured from the axis of the 
compressive stress S. 

If r>>a, S, and Tg approach the values given in 

equation 2-9. At the edge of the hole r=a we find 
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It can be seen from (2-14) that Sg, is greatest when @ = 7/2 
Br 3.8/2 (i.e. when cosz@=-1) that is at the ends m and n of 
the diameter perpendicular to the direction of the 
Ron. At -these jpoints (5,),5, = 3S which isthe 
maximum compressive stress and is three times the uniform 


stress S applied at the end of the plate. At the points p 


and gq, @% is equal to and 0 and we find Sg = -S. Se there 
is a tensional_ stress, resulting from the applied 
compression, in the tangential direction at these points. 


For the cross section of the plate through the centre of the 
hole and perpendicular to the x-axis, @ =7/2 and from 
Equation 2-13 

ta SeOeres | = S/2.(2ta*?/rt?+3a*/rt) -cmn erm neem nnn 2-15 
It 1s clear that the effect of the hole is of a very 
localised character, and aS r increases the stress 
approaches the value S very rapidly, within a few hole 
Giameters. The distribution of this stress S, as a function 
of r is shown in the Fig. 12 by the shaded area, regardless 
of the relative isotropy of a material. Having the solution 
equation 2-13 for compression or tension in one direction, 
the solution for compression or tension in two perpendicular 
Girections can be easily obtained by superposition, of 
stresses given by (2-13) with (@+90) for the angular 
co-ordinate. 

A closer approximation to the geological situation is 

shown in Fig. 13 where unequal compressions S and Ss, _ S?>s, 


are applied orthogonally to a plate with a small hole in it. 
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Figure 13.... A plate submitted to two compressions S and 5s 
with a small circular hole in it 
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At the hole boundary r=a, it follows from the superposition 
that 

See oe Od oon oe ons) COS2 OS 2-16 
and Sg is maximum, when cos2@=-1 with value (3S-s) when @ 
NJ ot 2tat msn), and (minimum when | cos2@ -=1. with 
VODMCmIS STS) ewhenay = 90) tee Cato, Cd). 

Where measurements have been made, horizontal principal 
Stresses in the upper crust range from equality S=s to about 
S=4s or more (McGarr and Gay, 1978). The values of the 
horizontal stresses across the principal planes in the 
neighbourhood of the borehole have been calculated by Gough 
and beblw’ GlOBZ)perory Variows welative evaiuesvtopetne S/s 

ratio and are shown in Figures 14 and 15. Figures 14 A and 
B show the individual stress distributions. The principle of 
the superposition of the two parts of the stress field is 
illustrated ine “Fig. 14 C forthe case in which'S/s =1.0 . 
For S alone the tangential stress at the walls of the hole 
varies from a minimum value of -S (tensile) across the plane 
parallel to the S-axis to a maximum of +3S, across the plane 
normal to the S-axis. With the superposition, the stress 
field has radial symmetry and the tangential stress at _ the 
wall of the hole is +2S as can also be deduced from the 
equation 2-16,,.From equation 2-16 when Sg = 25 = 2s, it is 
independent of @ and so failure may occur but without a 
preferred ~azamuth.Asmbhe! ratio Of) S) toe) si _approachesm fone, 
the determination of the magnitudes and directions of S and 


s becomes less reliable because of the lack of resolution. 
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Figure 14.... Superposition of stress states about a 


borehole due to two horizontal principal stresses of equal 


magnitude. 
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The resultant stress “fields for other ratios of S/s are 
Saown nuntet igure sebhiseeromeuthe icase when ::S/s i4= 33.0; the 
tangential stress at the walls of the hole varies froma 


minimum of 0 to a maximum of +8s (see Eqn. 2-16). 


ira 
if ty 


Si 15 ‘afodt: pike: & 4 


at- apa ae eee 


xi -” if 
ts rt mba x 39 


III. DATA COLLECTION AND ANALYSIS 
Data have been secured from wells distributed widely in the 
sedimentary basin of Alberta as shown in Figure 16. In 
analysing the records it was observed that in depth ranges 
characterised by oversized holes, the tool tended to rotate 
in some parts and ceased to rotate in other parts. On _ the 
other hand, the tool ceased to rotate in certain ranges in 

which the two orthogonal diameters were equal indicating a 

circular cross-section. Such observations led to _ the 

adoption of three criteria to be used in identifying 
breakouts. These are: 

1. The tool ceases to rotate, with the azimuth of the No. 1 
electrode approximately constant. 

2. The two orthogonal diameters are definitely unequal with 
the smaller diameter at bit size. 

3. The tool was rotating both below and above the breakout 
zone. This ensures that the rotation was stopped by the 
breakout and not by some other cause. 

The application of these three criteria causes the rejection 

of many elongated zones,but ensures that only definitely 

identified breakouts are used in the study. Table 3.1 lists 
the breakouts identified by means of these three criteria, 
from the available log records. 

Representation and Statistics of Angular Data 

We wish now to combine the various mean azimuths of the 
breakouts in order to obtain a representative breakout 


azimuth for each well. We may regard an angular observation 
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Figure 16.... Map of Alberta showing locations of wells 
analysed in this study 
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ae SepOMNenOmiasCirelesot “unit radius. A single “observation 
@ (0° $< @ < 360° ) measured in degrees is then a unit vector 
andthe data can be; “described as circular -data. If the 
vectouers mot directed live if «thevangles @. (0? is gs) 180°) 
and 180°+@ are not distinguished, the data can be described 
as axial or non-polar data. Ungrouped angular data can be 
represented in two ways; (a) by points on the circumference 
of a unit circle , the same weight being assigned to each 
observation, or (b) by; drawingsthe radiije® a unit icirele 
obtained by joining the origin to the observed points on the 
circumference, (figure 17). 

Alternatively, angular data can be grouped by dividing 
the range 0 to 360 degrees into a certain number of class 
intervals. The frequency is then the number of observations 
within each class. The data can then be represented on a 
histogram similar to that used on a line, by constructing a 
block whose area iS proportional to the frequency in that 
interval on the circumference of a unit circle. This diagram 
is called a circular histogram. This histogram can also be 
unrolled so that it sits on a straight line divided into the 
Glassmintenvals The point. Of Cub BiLOr sUnrOLl ings the Munit 
circle needs careful selection. For data that have a mode (a 
preferred direction) it is reasonable to use a cut such that 
the centre of the linear histogram approximately corresponds 
to this mode to avoid division of the modal peak between the 
ends of the histogram. This linear histogram is preferred to 


the circular histogram mainly because it iS easier to 


a 
' 


A = ee¢nbtid 28 


_ 


fe t m» 2 "CY © gations St; 7 aur 


: : “vie é t att. 2F .. ey "Oe 
7 a ate (ae ba é ruven te = a: a a s 


7 
ee aN 


eye bud - pa ine 4 


: 


; 

Pi 

= 

- 
e! 
~ hal. 
_ 
or 
— 


ate anae. edd ‘ vi caanad 


3 ir giid ‘ea9 9 
kit Sot 
. ‘ " - ro) 
i ee ; . \ Beeb Yeiopes’ (Win wesadie tA 
- ’ J an #30 
> ri t 
ats “ar pit SAF . cn 
: a 2iT .#a8s 
? ; 
Te 3-3 ary al o2 at batt 
his 26 shreze Tat 7xo. wt th 
~~ oan eh coy 4 - oote ¥ att i Bai os s ‘po! 
} beSiv ib sai> tdotetde a fe Bote oh) peas e neltoi bees 
, ~ el of ee a 
ome Litas fot jus Te Aieod ent alsyv ost au 


a) shon so aad tert {eb to: + iit ee ete twiarko abeait's 


1 


- ‘ - _ 
BRAGGdeSe 1 3CS - “les 7 ial i AO Woe fe BSt ut ‘i : aE i 
| ee - ; i oe b ule es 03 a ees 
" g9 ceevted “— i aap oat ae ere Te eres 


) 


— 
_ 


. : - oy ne 
ao rs 4 os ie Ai be: Pree ee wap 
me es 7 xe > \ 5 7 : J ot os : 
it ee nnn eee mareodaid taal Che smerpoie te 


Ea aoe oe 
al reieAe ai) By 


| ee eae 
ay Wiser se) exe 


Figure 
OP=1 


RES Ba 


Representations of the itn sample point @ 


50 


Pi 


Jy va cig. 2 irpome aids 


ty 


D1 


evaluate. 

Another natural representation of angular data is a 
Rose wdiagram. In) this “approach we construct a sector, 
corresponding to each class interval, with apex at the 
Origin and radius proportional to the class frequency (and 
arc substending the class interval). A disadvantage of this 
representation is that the area of each sector is 
proportional to N’?, where N is the frequency, so that the 
data appear much better grouped than they really are. 
However, the Rose diagram is widely used in geology. 
A Measure of Location 

Pethewerrnsts vinstance: ite: 2S) tempting. storeuse “<cehe 
conventional measures on the line mover tl. Velabiieibllicus 
distribution. Let us consider a case where we have a_ sample 
of size 2 and the observed angles are 1° and 359°. The 
arithmetic mean and the sample variance give unreasonable 
results. Intuitively, however, we infer the mean direction 
and the deviation about the mean to be in some sense 0° and 
1° respectively. A sensible answer, however, results if we 
select the zero direction as the y-axis in place of the 
X-axis.) This) will then reduce’ the data sto, 2699 and 271°: 
Hence, the usual linear measures depend heavily on _ the 
choice of the zero direction and will not be appropriate for 
Ginrewlarlaustrebutrons. Let. Xj 5.6. .... ,X,be n observations on 
Hemera tance etek. lees. sels cuisine pane Prepresenc ene .Same 
observations when the distances are measured from O' instead 
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Suitable measure of location on the line we should have 
L(X,',....,%,')= L(X,,....,%,) - a -oo crn = ao 
which implies that the position of the point whose 
x-coordinate is L{X,,...,%,) remains invariant under the 
choice of origin. Similarly it is desirable that the measure 
of circular location should not depend upon the choice of 
the zero direction. 
DeLee 7, De the. angles, obtained: from @,,..6.,0, ‘with 
respect to a new zero direction OA. If angle XOA = a, then 
we have 
ree on NO. 6 On 8). MOO ieN an a ao a 3-2 
The contrast between Equations 3-1 and 3-2 indicates the 
different natures of the problems on the line and on the 
Citcie., 
The Mean Direction 


suppose Bis the point on “the circumference of a unit 


exrcie ‘corresponding to the angle © , 1:= 1,...,n. Then the 
mean Givrection @ of @,,....,@,1S defined as the direction of 
the resultant of the unit vectors OP, ,...... HOR. 7 The 


Cartesian coordinates of P are (cos@. ,sin?, ) so that the 

centre of gravity of these points is (€ , S ) such that 
Clin 80s 6 Se 94/8 5i00, |). ---- Sac 
If R =\(€? + §?) --------------------------------------- 3-4 

then R = nR is the length of the resultant ane @ ise ‘the 

solution of the equations 

C = R.cos@ : S = R.sin@----------------- rr 3-5 


It can be shown that @ has some desirable properties as a 
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measure of location.(see Appendix A) 
The Circular Variance 

Bet, be the point corresponding to on a unit circle 
and a be a fixed direction. If we assume initially that «=0 
and suppose that P is the corresponding point on the circle 
then a measure of the circular dispersion between P and PR is 
the smaller of the two angles that OP makes with OP, say B. 


(see Fig. 17). We have 


Be aL a i ah ein map ne es Sie nea 3-6 
Cosh. AS van ancreasing function of £.. Let 

El GOSt eee yet pin Rena aera Hae | 
be a ymeasure of the dispersion of ~the ‘points  P. . After 


Shifting the zero direction to a using 
0." = (0..— a) modu27 , 
we can write 
Te ee es OG ana) ee ie ein ie 2-8 
The dispersion D is minimised at a=@ . If we equate the 
derivative of (3+8) with respect to a to zero we have 
Zsin(%, - «)=0, so that the dispersion is smallest about @ 
(see Appendix A Equations A-1 and A-6.). This is similar to 
the expression for the ordinary sample variance. D can be 
written about @ as S, where 
Sa edly emercesiie- Oa) G—e Cae gegeresn as ~~ eS 3-9 


which can-:be simplified to 


S, is called the circular variance and R is the mean 


resultant. vectorumagnitude, Frome (AcS) , eit \cangbé seen tthat 
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Se. ?S 2MVariant under a change of the zero direction. Let R 
be the length of the resultant of the vectors OPe Sites 7 R= 
Auk so lthat 

S, = 1 - R/n ----------------------------------------- 3-11 


Thus we see immediately that 


Since R tends to n for tightly grouped unit vectors, clearly 
S, tends to 0 for tight distributions and S, tends to 1.0 
for random distributions. From (3+12) S, takes values in 
(0,1) unlike o* for linear data whose range as. (0,infinaity):. 
An appropriate transformation (first suggested by Von Mises 
in 1918) of S, to the range (0,infinity) is given by 
g=4(-2.109, (1 - S,)) ---------------------------------- 3-13 
In (3-13) it has been assumed that the range of @ is (0,27). 
If the range of Y is ({0,27/L) then we define 
AE Ges Weys pas ta Din 5.2. | 
For grouped data we have 

Greedy nob cos@. Si ti) fs. Sings pi -vand oR Je (C7 + 
Seay 
Then cos@ =C/SR; 21nd =oS/iR.. See to - R 
To obtain @ it is convenient to use 

Ba! 38 -S > 0, C >.0 (ist quadrant) 


peut oC. 0. 7's Sx0 Lend quadrant) 


ISU 
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Ql 


+rif C <0, § < 0 (3rd quadrant) 


Q 
" 
Ql 


Gis OG * +9 if S.< 0 , € > 0 (4th quadrant) 
where @' = arctan(S /C ), and -7/2 < @' < 7/2. For axial 


(non-polar) data, in the range (0°-180°), we first double 
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the angles to maintain the periodicity. Thus we have @, '= 
2.0, and carry out the calculations using the equations 
already listed, and the frequencies f. corresponding to the 
mre-paints @ '. Suppose GC’; SS") @ " etc.-are«the quantities 
for @%,' corresponding to ¢,S , @ etc. then an appropriate 
measure of the circular mean @ for the OF can be taken as @ 
= $'/2 and an appropriate measure of the circular variance 
S. tor 2. is See ee In general if the range of @ is 
(0 , 27/L) we have 
See 51 ye 

The result of applying these statistics to the raw azimuth 
data measured at every two feet depth interval within an 
identified breakout is shown in the column headed MEAN in 
Table @.1 and also in columns 6 and 9 of “Table °3:;2. "Table 
3.2 shows the mean orientation of the various breakout zones 
in a well. 
Statistical Decision 

We are now near a position to make inferences on _ our 
data. Before making any inferences, however, we wish to look 
at the problem of the regression of one variable, in our 
case the azimuth, on another variable {independent 
variable), which is the depth in our situation. We shall 
then look at the correlation or the degree of relationship 
between the variables. We first consider the problem of how 
well a straight line explains the relationship between two 
variables. This requires the equation for the Least Square 


Regression Line. The least Square regression line of Y on X 
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TABLE Bie 


Statistics of breakout azimuths for wells listed in Table 3.1 
(All lengths are in feet except where indicated) 


WELL LOGGED ELONGATIONS 


LENGTH BREAKOUTS LENGTH N MEAN AZIMUTH(DEG) $.0(DEG) S.E(DEG) 
6-30-°46-17WS 3226-12106 80 7148 51 3080 1540 145.9 10. o.3 
7-9-49-24WS 1418-13939 144 9357 7 2336 1168 156.7 Ziti °.6 
12-19-53-8Wws 6350-7827 11 420 10 262 131 123).8 22 0.2 
b-45-A-94-P-14 5650-6089 13 103 6 32 16 60.8 Js 4.6 
6-12-50-11WS m 461-2874.5m 18 243 12 87 46 126.3 2. 0.4 
14-20-50-12WS 1230-9558 28 831 24 418 209 $22.4 4. 0.3 
4-26-S1-11WS 8484-9123 2 56 2 52 26 119.8 ES 0.7 
3-32-47-13WS 9566-11494 28 758 23 526 263 1345.5 ue °.5 
6-30-48-12WS 10011-11630 13 316 12 242 121 133.9 4. 0.4 
10°25-45-15WSs 10450-11335 13 298 1 204 102 148.3 35 0.4 
Bae as a 9000-13452 15 514 13 292 146 124.9 4. 0.4 
| 10-27-45-16Ws 13500-15407 20 2624 14 1342 671 135.3 ree ©.1 
13-9-S2-8W5 1526-7420 36 3098 20 1462 730 137.4 7. 0.3 
8-8N-51-39WS m 2200-2451.5m 1 78 8 30 16 134.7 3. (ekas 
6-20-50-15WS sseee 46 1330 4) 976 491 132.5 3. 0.4 
3-36-49-23WS 1507-18002 143 14518 84 3344 1672 151.2 10. 0.3 
5-9-S1-3WS 5550-8024 4 2455 2 340 170 167.8 25 0.2 
10-7-60-14WS 7500-8736 16 450 WW 340 170 154.3 8. ©.7 
2-29-59-22WS 7700-12864 72 2717 41 1880 340 135.0 10, 0.3 
14°21-37-11WS 5084-12040e 37 SO16 30 3060 1530 151.0 4. o.1 
7-7+86-1W6 795-7213 30 2483 4 190 95 114.3 3: 0.4 
9-18-83-6W5 5000-7675 12 876 10 710 355 115.4 4. 0.2 
5-5-35-11WS 14205-18000 33 1560 12 $40 470 170.7 8. 0.4 
| 15-29-115-3we 3500-4644 29 2215 12 690 345 °.6 10. °.5 
S-26-11S-4W6 4368-5205 _—s55 166 3 60 30 PSs seas 37 o 
12-24-82-6W6 4100-7073 19 1562 10 910 ass 105.8 6. 
11-29-82-SW6 4100-7311 31 1182 19 570 285 104.4 10. °.6 
3-21-112-6W6 4150-4802 6 146 2 40 20 159.0 °. 0.2 
| 13-35-112-10W6 4400-6200 22 948 12 360 180 150.8 ee ©.9 
{ 
3-1-119-9W6 4650-5774 18 2442 4 160 80 140.4 1S. a7, 
3-11-121-7W6 4800-5624 4 90 2 30 15 160.0 °. 2.0 
12-33-116-6W6 4250-SO78 38 2385 8 180 30 140.7 zo. 2.2 
14-36-116-6W6 4200-5300 4 146 3 80 40 374 43. 2a 
14-9-118-8W6 648-6091 27 2236 4 140 70 taal. 3 ie 1.4 
9-15-118-8W6 5042-6142 8 256 3 130 65 1285.7 °. o.4 
2-36-118-8W6 4600-6135 10 628 4 140 10 153.0 3. 4.2 
5-33-109-8Ww6e 692-6024 5 239 3 160 80 103.4 14, 1.6 
7-18RN-109-7WE 4800-5888 12 196 9 140 70 13.0 6. 0.8 
10-17RN-110-39W6 S400-6612 "1 1350 s 360 180 47,3 oe 0.9 
2-20-81-9WS 4000-5848 7 620 5 330 165 115.3 ©. o.1 
4-20-81-9WS 622-5855 28 1730 17 390 495 123.1 8. °.4 
5-3-81-9WS 3200-5871 20 1610 9 620 310 97.8 4. 0.3 
4-15-81-9WS 4900-5927 7 446 3 70 35 116.6 2s 0.4 
‘ 
| aS-5-82-9ws 5100-5590 4 98 3 40 20 67.1 s. Noe 
| 3-32-34-8ws 4000-9430 15 3034 9 2020 1010 133.3 T 0.2 
| §-19-37-9ws 7100-10263 41 1866 39 1730 865 121.9 2 O71 
A2-13-34-11WS UNSPECIFIED 9 182 7 128 64 138.3 Ss. °.7 
7-17-°33-7WS 794-8200 31 3460 17 2204 1102 136.3 13). 0.4 
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where a, and a, are obtained from the normal equations 

an. = Begun + (a,y. 2% 

ERY = a4. DK + ay. hk? eae ee eee mee sree ee ro Ball ites 
which yield 


(Sy)(dox2) - (Hx) (xy) 
N2x* - {)Ux}* 


ao 


Wek > aes. we have 9 positive ov .< ey, 
Nexen> oxi? 
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a, 1S known as the regression coefficient of Y on X. In our 
case Y represents the azimuth of a breakout and X the depth. 
To evaluate the significance of a regression coefficient one 
must consider the scatter or spread of points about the 
regression line of Y on &. This 1S given by the statistic 
Standard Error Of Estimate of Y on X.(Sy,x). Suppose we let 
Yess represent the values of Y for given values of X as 
estimated from (3.14), then 

Sy,x = {20 v- Yost) ?/ n} Gate e eg) 7a ety a ot ee ot ae eee et 
Equation 3-17 can be written as 

Dor eae ee a ae en Se ee Se S548 
This standard error of estimate has properties analogous to 
those of the standard deviation in the sense that if lines 
are constructed parallel to the regression line of Y on X at 
Fespective vertical distances Sy,“ , 2Sy,;x , and. 3Sy,x from 


it, it is found (if n is large enough) that there would be 


aay 
he Ett nt a nwt eH 
wAl ~~ - a 


> : ; 


0) i toupe nae eit. a bast 


sviissad « Xo nt ean eid atnesezge? 
| | . ; | | a = t se : ‘ 
j x » i IOs ves P se 3 i S, 2c oc : + . 3 A ‘imie oo - 52 on hes 
: ; 2 a. 


ee 
2 en 


aor 1 eetde. FoF iofsese ang vebignes: 


in snit Sulla 
8 


ray 


i+ vd evbe ai ail? sR a 


=+ 


soda kee te, Y pe a " Aone 


is. 2emdsy vio wot cat Joe ist saz aige'siqe’ i; 


| | Pe oe kent) ‘mont vere 
aes onal Bae ne = ens pieninilg |! Hit 2X sas 
| | ap gdh. os tae Thnk nek, 


5 7 | . - . rok 
aT Bleep mat nena (ree a ie omplaapeaael 
n | ons 
a eo eee ae ablinege’ Ban aagninty 
et y > v - ny i, 


a Sedna ada fi, Aagtebe 


ath anit otatas 


As orotic 


2 


included between these lines about 68% , 95% , and 99.7% of 
the sample points respectively. Now that a least square 
regression line has been fitted to the points,we would like 
tovknow (how well ‘a straight line fitsr the ‘points “or 
describes the relationship between the variables. In other 
words we want to determine the goodness of fit of the line 
fitted to the data. The quantity which does this is the 
coefficient: of correlation (4). ro ois a dimensionless 
Guantity which varies between -1 and +1. If ¥Y tends to 
increase aS X increases we have a positive or direct 
correlation. If Y tends to decrease as X increases we have 
negative or inverse correlation. If r 1s near zero it means 
there is almost no linear correlation between the variables. 
r can be expressed in several ways but for computational 


purposes we write 


\ inox2- (Sox HNDv2-(Sov) 


It is worth noting that a high correlation coefficient (i.e. 
near 1 or =i) does “not: “necessarily indicate a direct 
dependence of the variables. 
Tests of Hypothesis and Significance 

In order for any test of hypothesis or rules of 
decision to be good they must be designed so as to minimise 
errors of decision. Obviously this iS not simple at all 


Since, for a given sample an attempt to decrease one type of 
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error is accompanied by an increase in other types of error. 
However, one type of error may be more serious than another 
which may (call for ayscompromisesin favour of a limitation of 
the more serious error. The only way to reduce all types of 
errors 1S to increase the sample size which may or may not 
be possible. 
Level of Significance 

Rejection of a hypothesis when it should be accepted is 
defined as a Type I error. Acceptance of a hypothesis which 
Should be rejected is defined as a Type Il error. In testing 
a given hypothesis, the maximum probability with which we 
wish to risk a Type I error is called the Level of 
Significance of the test. It represents the probability of 
our being wrong in rejecting the hypothesis. 
The regression equation Y =a, + a,.X was obtained on the 
basis of sample data. We are often interested in the 
corresponding regression equation for the population from 
which the sample was drawn. 
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Pop pop 


If we wish to test the hypothesis that the regression 
coefficient of the population A, 1s zero at some chosen 
Significance, we use the fact that the statistic t has 
Student's distribution, with (n-2) degrees of freedom, where 

f= (ta, -Ayiesx): fee ere - (a,-Ay) ff(n-2)/ (1-r?)t 3-20 
Here Sy,x is the standard error of estimate of Y on X, Sx is 
the stendara deviation. Of X,. vand- ri is the ‘correlation 
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intervals for population regression coefficients from sample 
values. 

In applying @Equawions--_oal4.athroughasand including 
Equation 3-20 to our data we can then compare our t_ values 
to the Student's t-distribution table (given in many 
Statistical text books), with the correct number of degrees 
of freedom, and find at what level of significance the 
hypothesis is to be accepted or rejected. 

Table 3.3 presents regression results for 21 wells in which 
breakouts were meaSured over a depth range greater than 2000 
feet. Column 3 shows values for the regression coefficient 
a, for the sample and column 4 gives the standard error of 
estimate of the azimuths on depths. Column 6 shows the 
various t-distribution values and column 7 shows the degrees 
of freedom which can be used to find the level of 
Significance at which the test can be accepted by reading 
the appropriate t-value and degree of freedom from tables in 
text “pooks. In ‘column 5 we have the -values of the 
correlation coefficient. 

Results 

It is clear that the regression coefficient a, 1s inall 
wells small enough to have arisen with more than 5. percent 
pLObaDI ity, ein am population having a,=0..0f ‘the, 21.welis 
tpeteduan Taple 3.5 4 10 ehave a> Ovand. 1 i have.a ys <0, 1. We 
cannot reject the hypothesis that the azimuths’ of the 
breakouts are unrelated to depth and the variation observed 


with depth is only by chance 
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TABLE vd2a 


Regression of breakout azimuths on depths showing the 
regression coefficients a, (deg/ft) , the correlation 
coefficient r, the standard error of estimate Sy,x(deg), 
the t-distribution, and the degrees of freedom df. 

(Only wells with the entire breakout intervals extending 
through 2000 feet depth intervals are listed here). 


Well a, (deg) a,(d/ft) Sy,x(deg) t df 
6-30-46-17W5 Looe o 0 e002S2) 28737 0.183. 4-0. 059 619 
7-9-49-24W5 Hib. 9) 0.00112 45.7 0.087 0.038 tt24 
14-20-50-12W5 129-78 -0.00080 10.8 -05.052 70071 199 
10-27-45- 16W5 tyOgle -0.002390 4.1 -0.460 =0.085 680 
13-9-52-6W5 143.37 -0.00242 12.3 -0.203 -0.050 405 
6-20-50-15W5 ZOAatt -0.01416 18.1 -0.475 -0.345 459 
3-36-49-23W5 143.06 0.00042 28.7 0.062 0.018 iMhew 
5-3 =81-=9W5 95.82 0.00046 1029 0.024 0.008 308 
4-20-81-9W5 Weeese 0.00055 16.9 0.050 0.012 493 
3- 1- TAS ZOWS 17269 -0.01128 30.8 -0.489 -0.114 78 
12-33-116-6W6 144.05 =0 00577 41.2 = 0.269 -0.066 88 
11=29-82-5W6 35.66 0.01150 20.2 0.356 0.200 263 
i pa ah hl Be a 9) 15370 0.00856 ZU% 3 0.304 0.166 343 
9305.30" sW5 = 289.62 0.02845 1a 0.716 0.881 468 
9-18-83-6w6 16452)2 -0.00786> 1076. *=0'2306 ~-07154 348 
14-24-37 =4,0W5 eh BES -0.00011 bial “0037 -0.005 1838 
229750 722WS ie. 99 0. 00t7Ban22:9 0.115 0.055 938 
3-32-34-8W5 167.45 -0.00442 biccd -0.446 oO2 158 1022 
6-17-37 -SW5 130574 -0.00105 oxo Uw MS) -0.031 863 
R2-V3-34"1 1W5 65.50 0.00553 14.5 0.194 0.044 62 
i 1=33=7W5 72.14 0.01084 3iio 0.933 0.420 1073 
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IV. BREAKOUTS AND STRESS ORIENTATIONS IN THE WESTERN 
CANADIAN SEDIMENTARY BASIN 
Relation of Breakouts to Location and Lithology 

In Chapter III] the statistical distribution of breakout 
azimuths has been set out. In most holes, azimuths are 
closely grouped together, as shown by the mean orientations 
of the various breakout zones within each hole (Table 3.1), 
and the standard deviations listed in Table 3.2. In the 
holes logged through a sufficient range of depth to allow 
investigation of regression of azimuth on depth, such 
regression is not significantly different from zero (see 
Table 3.3.). Figure 18 is a map showing how the azimuths 
enange with | locations, Thes azimuth—-referred to is ‘the 
representative azimuth of breakout for the well as a whole 
obtained by combining the various breakout azimuths in the 
well given in Table 3.2. The central dot shows the position 
of the hole and the mean orientation of the breakout is 
Shown by the line through it. The number attached to the 
line is the mean orientation of the various sections in the 
well. 

We now consider the relation of the breakout azimuth to 
the geological age (i.e. period), formation and lithology. 
These data are shown in Table 3.1 under the columns 
Rormation. (Grn)  Geashi thology “(hithy, and Period (Age)), ‘and 
have been obtained from the sample descriptions for the 
wells. Sections under these columns which are blank indicate 


that the information was not available for those wells. It 
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1s evident, however, from these columns that breakouts have 
occurred throughout the sedimentary basins in rocks ranging 
from Devonian to Cretaceous in age. In some cases breakouts 
have been observed at different depth intervals running 
through the entire depth of the hole. 

A look at the columns headed formation and _ lithology 
revealS many interesting points. Breakouts have occured in 
formations like Elk Point, Wabamun etc. in the Devonian 
right to Belly River and Bearpaw formations in the Upper 
Cretaceous period. Within a given formation some sections 
show breakouts but other sections do not, for reasons which 
are obscure. BreakoutsS may occur in any position such as the 
top, middle or bottom sections of the formation. The pattern 
does not -appeane ho, tThemiregulatedrebys theevtypeceofhf the 
lithology. A word of caution must be sounded here. The 
sample descriptions from which our stratigraphic information 
was obtained merely stated the predominant lithologies found 
and only stated the traces present but not the percentages 
of these traces. A knowledge of the percentages of such 
traces would be useful in view of the fact that there are 
several diagenetic processes involved in the transformation 
of unconsolidated sediments to consolidated sedimentary 
rocks. However, there are formations in which only one kind 
of lithology, without any traces, runs through the entire 
depth range in which breakouts occur in certain depth 


intervals separated by uncaved intervals. 
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Although the elongated sections which did not satisfy 
our three criteria for selecting a breakout are not included 
in Table 3.1, a careful study of the ages, formations, and 
lithologies of these zones was carried out. The entire study 
shows that breakouts are not lithologically controlled. They 
have occurred through most parts of the Western Canadian 
sedimentary baSin cutting through limestone, dolomite, 
Shale, sandstone, anhydrites, cherty shales, siltstone, and 
varying combinations of these sedimentary rocks. At the same 
time, the intervening zones without breakouts have exactly 
the same types of lithology as already listed. It can 
therefore be safely concluded that the breakouts and/or 
elongated zones. are not’ significantly related to the 
lithology and age of formation. There is no evidence of any 
Significant relation between breakout azimuths and _ the 
properties of the rock. 

This agrees with the observation of Babcock (1978) ona 
smaller sample of breakout data. The logs used in this study 
give no indication of the dip angles of the various’ strata. 
This study therefore cannot add to the evidence given by Cox 
(1970) and Babcock (1978) that the breakouts and their 
azimuths are unrelated to the dip. 

There are two wells included in this study in which the 
two breakout zones selected using the three criteria have 
azimuths approximately orthogonal to each other. These wells 
are included in Table 3.1. Some other wells have one or two 


intervals with breakout azimuths differing by as much as 90° 
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from the dominant azimuths. Such wells have not. been 
referred to as bimodal wells. This study has shown that 
whatever their cause, a predominant direction of breakout is 
almost always present. Variations in this direction with 
position may represent local variations in the orientation 
of the horizontal stresses. 

Pote ewatera boeostre 7h. ac tne pressure oP) i ihe 
Grilling mud which fills the hole were also disregarded in 
the discussion in Chapter II. These have no effect on _ the 
crientations of the fractures (Gough and _ Bell,1982) 
nevertheless, they modify the magnitude of the stress 


differences. At points P and Q (p, q) of figures 12 and 13 
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Pee ca hace lal a kee ie = 0. 


Thus the stress difference becomes 


ThesMonr. circles fig.) 10 contracts in diameter by. BP» and 
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Loup The vchear “sarengehyt, Jrequizedkto, resist failure wild 
be reduced somewhat as a result of the shift of the Mohr's 
circle to the left and the nature of the envelope. The shift 
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stresses, in a biaxial stress field, needed to cause 


failure. 
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Hydraulic Fracturing 

In-Situ measurement of stress by hydraulic fracturing 
can both provide the orientation and magnitude of the 
Principal Componencs. -of “stress,” it ~S,)°is not vertical 
(Kehle, 1964 ; Fairhurst, 1964), provided the fracture 
pressure, the instantaneous shut-in pressure (ISIP), the 
pore pressure and the tensile strength of the rock are 
measured. Hydraulic fracturing is a petroleum engineering 
technique carried out at depths of several thousand feet to 
stimulate production of oil bearing horizons. An 
understanding of the basic principles which govern the 
development of the fractures iS a prerequisite to its 
successful application. It is assumed that the hole is 
drilled parallel to one regional principal stress usually 
the vertical stress S,, which 1s assumed to be equal to the 
lithostatic load p.g.h. The stresses S and s normal to the 
hole axis are then the other two principal stresses and are 
horizontal. The pre-existing regional stresses produce 
stress concentration close to the borehole with the maximum 
value of (3S-s) and a minimum of (3s-S) -in ‘the tangential 
direction. The smallest tangential stress of (3s-S) at the 
ends of the diameter parallel to S may be greater than, or 
equal to, or less than zero depending on the S/s-ratio. If a 
section of the hole is sealed and pressurized to a pressure 
p, a second stress system will be superimposed on that just 
described. In the region ciose to the packers, Kehle (1964) 


assumed that the packers are held in place by a band cf 
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uniform shear stresses and made an analysis to determine: 
these stresses. His resultS indicate that high tensile 
vertical stress of the order of the applied pressure (-p) is 
induced in the immediate vicinity of the packers, and in 
about 80% of the central region of the interval a tangential 
tensile stress is generated. If however inflatable packers 
are used to seal the hole a normal pressure will be exerted 
resulting in a considerable reduction of the axial tension 
and will require a much higher pressure for the initiation 
of horizontal fractures. Before such high pressures are 
reached a fracture will be initiated vertically (Zoback et 
al, 1977). As the pressure in the hole is increased a 
fracture will be initiated at the wall of the borehole when 
the tension generated by the fluid is sufficient to overcome 
the combined effect of the regional compression (3s-S) and 
the rock strength. The plane along which a fracture will 
commence will be that across which the resultant of the 
effective compressive and tensile strengths are first 
reduced to zero. In the case of a smooth cylindrical 
wellbore this plane must be vertical and perpendicular to 
the least principal regional stress. If the vertical stress 
is the minimum principal stress a horizontal tensile 
fracture might be expected to form (Kehle,1964; Hubbert and 
Willis, 1957) However, Zoback et al (1977) have argued 
against this. This contention of Zoback et al 1S supported 
by laboratory experiments reported by Haimson and Fairhurst 


(1970). Haimson (1976b) suggested that if S§S = §, then 
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induced fractures initiate in a vertical plane and then 
become horizontal as they propagate a few hole radii from 
the wall into the undisturbed regional stress field. 

The breakdown pressure FR is recorded. If the pump is 
shut off immediately and the hydraulic circuit is kept 
closed an instantaneous shut-in pressure, a pressure that is 
just sufficient to hold the fracture open, is recorded. The 
theory derived by Hubbert and Willis (1957) and Kehle (1964) 
for fracture around a pressurized borehole is used to relate 
BP the breakdown pressure, the instantaneous shut-in 
pressure ISIP, to the in-situ principal stresses and the 
tensile strength, T, of the rock. For a vertical borehole, 
the tensile fracture should be oriented ina direction 
perpendicular to s, the least horizontal principal stress, 
at least within the immediate vicinity of the borehole, and 
the magnitude of s is equal to ISIP. S, the maximum 
horizontal principal stress, is determined from the equation 

ae pce eg eI NI a a aie ace el ara ee od din ce ch So 
where p is the static pore pressure in the rock surrounding 
the borehole. T can be meaSured in the laboratory and p in 
the field. From the discussion so far it 1s clear that the 
hydraulic fracture will generally take place parallel to the 
maximum horizontal principal stress. Breakouts on the other 
hand occur parallel to the lesser horizontal principal 
Stress or perpendicular to the larger horizontal principal 
stress. It is therefore, not surprising that fractures 


detected by the 4-arm dipmeter tool may occasionally show 
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two trends one being a major one, if mud_ pressure _ has 
occasionally produced inadvertent hydraulic fracturing. 

The question to be asked will therefore be, how the 
hydraulic knactvuring comes inoogiplaye. Pnoior ait of fehe 
deliberate application of hydraulic fracturing by man, its 
occurence has been detected during water flooding and 
squeeze cementation processes. Rocks have been accidentally 
fractured during grouting treatments and when high specific 
gravity muds have been used in drilling wells. Furthermore, 
it has been observed (Von Schonfeldt and Fairhurst, 1969) 
that formations without any measurable permeability under 
habuncalh, conditwons, Have).a noticeable effect of rate of 
loading. At a low rate of application of pressure, the 
formation breakdown is lower than at a higher loading rate. 
Such a loading rate may be likened to the varying rate at 
which drihbing is done depending on the commercial 
importance of the zone being drilled. Accidental hydraulic 
fracturingadduring isthe 18drilbinghastagepowrthdhighaspeerinc 
Gvuavney mudicanmtheneforernottbengruled outiawibticns wonth 
mentioning also that although the vertical stress 
approximates the overburden, thus implying a reduced effect 
of accidental hydraulic fracturing, we should not lose sight 
of the fact that the oil wells are restricted to sedimentary 
layers. Sediments are sheet-like deposits of various 
competency stacked one upon another. Layered sequences of 
unlike materials transfer stress much differently than 


homogeneous and isotropic materials, a phenomenon known as 
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arching. Arching may also involve the reaction with time 
between beds of different competency and hence disrupts the 
uniform stress distribution with depth which is assumed 
under homogeneous and isotropic conditions. Such an effect 
creates regions of both high and low stress concentrations. 
As a result of the alteration of the stress distribution by 
anchimotwe camgexpeatdthatdinethe! sdmekteibs field, inforcva 
Given horizon breakdown pressure can vary significantly as 
evidenced by the results obtained in the Gulf Coast area 
where vertical fractures have occured at injection pressures 
far less than the total overburden pressure ( Scotts, 
Bearden, and Howard, 1953). 

Another point which may need examination is the 
assumption that one principal stress is vertical and that 
the borehole is drilled parallel to this regional principal 
stress. The assumption that the hydraulic fracture 
propagates perpendicular to the least principal stress ( 
Hubbert and Willis, 1957) has been supported by agreement 
between hydraulic fracturing p geological and 
seismologically determined stress field indicators (Zoback 
and Zoback, 1980 ) However, McGarr and Gay (1978) found that 
in South Africa the directions of the most nearly vertical 
Pummeipal stress italle within = as icincle. tof, 30° about: the 
vertical. They also mentioned that stress meaSurements in 
deep mines in Canada, the United States and in Australia 
Support significant departures from the assumption that one 


of the principal stress directions is vertical. This implies 
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that the two initially assumed horizontal stresses are not 
Eluly horizontadakRecognitiom of such a.fact will also help 
in accounting for the irregular nature of fractures. 
Breakout and Stress Orientations in Western Canada 

Figure 19 shows the directions of the long axes of 
breakouts (or the dumecenon, of the Wleaste principal 
horizontal stress) and figure 20 indicates the inferred 
Orientations of the maximum horizontal principal stresses 
for the results presented by Babcock (1978), by Gough and 
Bell (1981) and in this study, for the West Canadian 
Sedimentary BaSin. It can be seen from this map that 
breakoutS can serve as a useful method for mapping the 
various tectonic stress zones based on the stress 
trajectories in any given locality where directional 
spalling of the walls of a well are observed. The absence of 
Significant variation of azimuths with depth, wherever these 
phenomena are observed, may imply that such large horizontal 
principal stresses are not necessarily limited to the top 
crustal sedimentary rocks but may be continued down into the 
deep and older Precambrian igneous’ rocks of the earth's 
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EXPLANATION 
Mean azimuth of 
5 or more breakouts 
No. of breakouts 
Mean azimuth of 
4 or less breakouts 
No. of breakouts 


Data from Babcock (1978) 
and unpublished 


Data from Gough and 
Bell ( 1981) 


Mean breakout azimuth in this 
study and no. of breakouts 
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Figure 19.... Map of Alberta showing the orientations of 
the smaller principal horizontal stresses as inferred from 
the breakouts 
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EXPLANATION 


Maximum horizontal 
stress orientation 
inferred from breakouts 


Minimum horizontal 
stress orientation 
inferred from breakouts 


Breakout data from Babcock (1978) 
and unpublished 


Breakout data from Gough 
and Bell (1981) 
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0° 


Inferred directions of the larger horizontal 


principal stresses for the results presented by Babcock 
(1978), Gough and Beti (1981) and in this study 
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V. DISCUSSION 

The Physical Significance of the Breakouts 

Babcock (1978) concluded on the basis of his joint-sets 
Study in Alberta (1973) that breakouts are a result of the 
drill encountering steeply dipping zones of pre-existing 
faults "or ~-jJoint=sets) ‘which’ may® lor: may »enot «che. :open: 
Preferential breakage, in his view, is then controlled by a 
near-vertical joint which intersects the well and which in 
carbonate rocks may be solution-widened. In his view, the 
parallelism between azimuths of elongation and regional 
joint sets in Alberta and the consistency of the azimuths of 
elongation within, and between wells in varying lithologies 
all indicate an origin caused by the drill encountering a 
PracturesorSzonecol fractures 

However, his study of joint-sets revealed four 
concentrations of joint directions on the surface, (NW, NE, 
N, and E). On a non-polar or axial data basis these four 
joint-sets directions may be represented as NW-SE, NE-SW, 
N-S, and E-W. Babcock's study of breakouts showed only one 
Significant concentration of subsurface breakout direction 
that is NW or NW-SE. Bell and Gough (1979) argued that on 
Babcock's hypothesis the other joint directions would be 
expected to be represented in the breakout azimuths,but they 
are not. In addition Bell and Gough doubted that the joint 
directions at the present surface would necessarily be 
parallel to those at depths in excess of 2 km at which 


breakouts occurred,.in view of the fact that sediments 
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involved in his study might have been laid down in varying 
tectonic settings between Devonian and Cretaceous times and 
that joints in general arise in several ways; such as _ soon 
afteretdeposituon Sof san sedimentaryarrock, or during’ later 
deformation, jouvas iteundergoesmeerosion: ifheses facts) “led 
Bell and Gough to put forward an alternative hypothesis to 
better explain the cause of breakouts and, perhaps, also the 
cause of surface joints of Babcock's System 1 (NE and NW 
Strikes). They therefore proposed that in a general stress 
field, with large and unequal horizontal stresses, the holes 
themselves concentrate the stresseS SO as _ to produce 
Subsurface breakouts with the breakout azimuth parallel to 
the smaller horizontal stress and normal to the larger 
horizontal stress. The analysis of such stress concentration 
by a circular hole was first found by Kirsch and has_ been 
discussed in chapter II]. 

From the present study, it is observed that the 
Significant breakout azimuth or the preferred breakout 
Orientation is SE, that is NW-SE when expressed axially. 
ThusT Conly  oned jot rt thes fourgwsurface] lLjomntrdinectionsi as 
Significant, although a few breakout zones tend to show the 
NE or NE-SW strike. In short only the NW-SE set of Babcock's 
System 1 of surface joints is observed as breakouts. On Bell 
and Gough's hypothesis that breakouts are stress controlled, 
the observed irregularities have been fully accounted for. 
Babcock's hypothesis fails to explain why breakout azimuths 


do not exhibit the other surface joint directions. Babcock's 
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study also revealed that breakouts are unrelated to depth, 
letbhologye and? thesrages ofgntheortormationst it asmiurther 
observed by Babcock and in this study that breakouts are 
discrete, beginning and ending Sharply, independently of 
deprn,siatholooyaandeagervor formation. Iltedsivery ditfiicult 
to understand how breakouts can begin and end sharply, 
within and between formations, separated by Similar 
lithologies with no regard to depth. In view of the various 
diagenetic processes involved in the transformation of the 
deposited sediments to consolidated sedimentary rocks from 
the Devonian to the Cretaceous age, which is more than _ 10° 
years, one wonders how surface joints can be a reflection of 
subsuriace @90Oints:s InenfiactSufactors like cementation, 
compaction, recrystallisation, and many other processes 
involved in diagenesis coupled with the high temperatures 
and duration of the tectonic stresses will modify the 
response of the rock to stress, while any failure remains 
brittle. Suchrwaraables! (may t lead’ oto “daillures, inu>some 
intervals in the sedimentary column, whereas the same 
material in adjoining area may not fail under the existing 
stress field. Consequently Bell and Gough's hypothesis can 
account for this discrete nature of breakout occurrence in 
wells. The presence of pre-existing fractures would also 
result in enormous circulation loss, which is however very 
small in this region. It is known that extensive circulation 
loss often leads to formation damage. Such damage often 


requires that hydraulic fracturing be carried out, this time 


peta ae ey ve Sears | 
ss whantie:s adh All morn cod a tov 
to ql ainebieegabat velerradia pets bie nee oe Bie 
itontiae emeae ad snot emiet ae eonteds) 

sigqiuad® She Sts mead, ‘hi ee ae 
ietinke” ae betatiqes ees sat on s.. at 
auoisev sda 40 qate mF ‘Wha os SiRget onthe: @ ; 2 gil” a ' 
aaa Ko Abi dgnesensts ane gh bertewnt: lanpesso1q: ‘oitonspai® + 
mond Py yisineminee satsbiloanns ot wtih” bes kaogeb a 
‘Ot; nets +ea0m'es dix beh ‘vans eugenssea? sats, ‘os astnoved: ea! | . 
| to ‘agbioatier Sad, 16D: pac ellbiaraeh weet aaebaow dno ame a a? 
aeiserneiies  o0b2)- siosany Sel = OL, L, | want, soetwuedue Jee 
apekerowy tego nin | itp. noi geni Loatayasea 0433808 : i 
aetusstecmes. ap dts ots ujiu balques ajastiegeit ai, boviovnt: rs 
sid yithou : ch esetasse aio soe silt. Do! ot semub 9008 
eciames sauhigl. gts, sitdv: egetse oe cers ent i eenogest a ey 
sme Ai Sy hSSe. eek yam! aeteilag polite | eigatad” fe 
‘aaa en? euatiedy sii tmenetbee: ais awh ahevietad 4 
oir esain ais ena £3 2 por, jem 9ts, ealNeDtb® i ete 
10. giasiioggt aniiguod m6 ers seins. stoid dat — ne 
fit eats 22300 Rodag ss te su sa8 3 alae 98 © ni 
| oes bTugk aviuitnent iui 
gna gewewortied: ete oa 
niatamtunsie: av beriene, a. 


nw 
= 


AY. ie ay ns i " a ad i 


92 


not just aS a means of providing an increased permeability 
path into the otherwise tight formation, but as a means of 
creating a conducting channel through the damaged zones into 
the virgin formation. The absence of such effects as high 
circulation loss, damaged formation resulting from the 
Grilling fluid escaping into the formation, and the complete 
absence of massive hydraulic fracturing processes in Alberta 
only go to support the point that pre-existing fractures of 
size as implied by Babcock are unlikely in the subsurface of 
Alberta. The only hydraulic fracturing process reported is 
the one at West Pembina (MacLeod, 1977). The fracturing 
direction 1S in agreement with the known theory on hydraulic 
Bractumiag andeithusonkendss supportudeo Bell thandprGouginss 
hypothesis (Gough & Bell, 1981). It is therefore, the 
present author's opinion that in the absence of positive 
extrapolation of surface joints to subsurface ones, 
Babcock's hypothesis does not account for the azimuthal 
alignment of breakouts. It is therefore concluded, subject 
to comparison with other stress meaSurement techniques in 
future, that breakouts are a result of the stress 
concentration by the borehole, in a general stress field 
having large and unequal horizontal stresses, at the walls 
of the hole resulting in the failure of some of the rocks. 
This is the hypothesis postulated by Bell and Gough. 

Based on the Bell-Gough hypothesis, we find that the 
greatest horizontal principal stress trends N 25° E to N 60° 


E; with othe majority occuring around N 42°°E, and. that» the 
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least principal horizontal stress is directed N 30° W to N 
65° W, with the majority around N 48° W throughout most 
parts of Alberta. Sbar and Sykes (1973) have reported the 
greatest principal horizontal compressive stress to. be 
teendingskast Co*‘Northeasttan BastermeNonth Americainand) jon 
the intermountain seismic belt in Montana, Idaho and Utah 
the minimum compressive stress trends East to Southeast. Our 
result thus fits closely to the general pattern of principal 
Stress orientations reported in the older craton of the 
North American continent. 
Implications of the Study 

Measurement of stress by in-situ techniques can provide 
both the orientations and magnitudes of the principal 
components of stress. Such measurement of the stress field 
within the crust can provide, perhaps, the most useful 
information -concerning the forces responsible for various 
tectonic processes, Such as earthquakes. For example Sykes 
et al (1972) reported the occurrence of a series of small 
earthquakes generated in 1971 by the high pressure injection 
of fluids ina salt recovery well near Dale, about 10km east 
of Attica in Western New York State. When the high pressure 
injection ceased seismic activity dropped in three days from 
a rate of about 100 events per day to a rate of a few events 
per month. The triggering of the earthquakes was attributed 
to the increase of fluid pressure in the rocks under high 
tectonic stress. Healy et al (1968) also reported the Denver 


earthquakes. The Rocky Mountain Arsenal disposal well, which 
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was drilled through 3761 metres of sedimentary rocks of the 
Denver BasSin near the city of Denver, Colorado, was intended 
forathetdirsposalGefawastevilurdsh The’ Ginjection of flurd 
which was started from March 1962 to September 1963 was at a 
bavesortabouth2).. ve10e Mites’ per Cnonth. PSPLuid® injection 
ceased from October 1963 to August 1964, but was resumed 
under gravity from September 1964 until March 1965 at a rate 
op about Meessnet 10°a Ireres” “per ‘month! *The raté “was *then 
increased to about 1.7 * 107 litres per month until February 
1966 when the exercise was terminated because of growing 
evidence of itS connection with local earthquakes. Evans 
(1966) used the epicentral locations by Wang (1965) of 
earthquakes near the well and showed a correlation between 
rate of fluid injection and the earthquake frequency, with a 
lag of the earthquakes in the order of a few weeks behind 
the injection. The events which initially had epicentres 
less than 6 km from the well, were observed to continue with 
migration of the foci to the northwest. Healy et al examined 
the probability of a chance association, in both time and 
place, between the fluid injection and the earthquake swarm 
on the basis of the seismicity of the region and estimated 
thisiteprobabi lie seto be above. ian 255 irri bien fe Healy set al 
(1968) therefore explained both the occurence of the 
northward migration and the occurrence of the three large 
earthquakes, among a series of earthquakes, of magnitude 
>5.0 to the outward diffusion from the well of a pressure 


front in the pore water. This diffusion would imply that the 
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volume of rock affected by the increased water pressure 
likewise grows and this increased volume is thus associated 
with the the large earthquakes in a large volume of major 
StraingssBut itecrrs) well) knowns that ‘athe increase. of the 
pressure of the water in pores and cracks produces no- shear 
stress and therefore cannot cause an earthquake. The 
induction process of the observed earthquakes was attributed 
to the presence of a nearly critical initial stress which 
could eventually cause failure of rock by triggering through 
by increase of pore water pressure as the Mohr circle is 
shifted to the left. This explanation is based on the 
Mohr-Coulomb failure theory, and the Hubbert and Rubey 
hypothesis of effective stress (1959) ( Healy et al, 1968; 
Gough and Gough, 1976; Gough, 1978). 

The Rangely Experiment carried out by the United States 
Geological Survey at the Rangely Oilfield in Northwestern 
Colorado (Raleigh et al, 1972, 1976) was designed to verify 
quantitatively the explanation of the induced seismicity at 
Denver. There was an excellent agreement with the 
hypothesis. Raleigh et al (1976) also pointed out that this 
verification might have implications for the control of 
natural earthquakes at least in cases where it might be 
possible to pump water in and out of an active fault. This 
could possibly relieve the stress by first reducing the 
water pressure in wells at two points to strengthen the 
fanltipsat athosetlipoints7s andietheniafinallys triggerang*an 


earthquake between the two points by injecting water there. 
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Knowledge of the state of stress 1s also critical to 
the design of underground excavations for mining and for 
nuclear waste disposal ( Jaeger and Cook, 1969; Gough, 1978). 
The massive hydraulic fracturing of formations in oil and 
Gasitueldsstoestumulate production rise ‘another: eiieldy for 
which knowledge of the stress field at depth 1s very 
important. Application of stress measurements to the 
solution of problems in tectonics is not as (straightforward 
as in engineering problems. The engineer iS concerned with 
the stress field affecting the rock, whereas’ the 
geophysicist or geologist attempts to deduce the processes 
that might have caused the stresses. Thus the importance of 
the knowledge of principal stress orientations cannot be 
overemphasized. The use of breakouts to add principal stress 
Orientations to the existing information on seismicity and 
tectonics will significantly increase the data base relevant 
to intra-plate earthquakes. Maps such as that shown in 
Brqureer20 sof sithearcorientatironnei the homizomtalh principal 
stresses, will be of great value to engineering purposes and 
in aqcaueion swrhll thavesvimportantiocimplications: efor) ithe 
geophysics and structural geology of the area. 

It must be emphasized that any attempt to relate the 
Stress field to tectonic processes, must be based on a 
correct interpretation of the stress field. A variety of 
postglacial geologic features such as folds? faults? 
pop-ups, and other rock squeeze indicators may show high 


compressive stress in a region. Residual stresses may exist 
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timed nockhpwaceoraing (Lomi tseghistonyve0die processes » such { zas 
bunpal, @urthification, @edenudation, shéeating, “cooling, Tand 
past tectonic events and may persist to some extent after 
the rock is freed of boundary loads ( Friedman, 1972). Such 
residual stresses may complicate the interpretation of 
Srresselobservat tons withoweversnnettewould bhe reduiircullt to 
attribute a consistently oriented stress field over a_ large 
area, aS in the entire Western Canadian sedimentary basin, 
to residual stresses. The observed stress field almost 
Gertainly ‘seflects tractions nowractingsonrthis part (of »the 
North American plate. The recent review by Zoback and Zoback 
(1980) shows that northeast-southwest horizontal larger 
compressions dominate the older craton of North America, as 
in the Northeastern United States ( Sbar and Sykes, 1973) 
and the Western Canadian data fit well in this overall 
picture. 

A knowledge of the orientation of the horizontal 
Stresses 1S of great importance in the planning of a 
hydraulic fracture program. Near-surface ore bodieS can _ be 
economically mined by open pit methods. Below the practical 
depth of open pit mining, new mining methods must be sought 
and developed to make deep ore bodies available to economic 
extraction. Solution mining has become a major method of 
extracting subsurface evaporite minerals because it is both 
economical and efficient. Hydraulic fracturing, apart from 
the fact that it is the only method currently in use that 


enables measurement of stress to be made at large distances 
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from a free surface, 15S used to develop permeable surface 
areas exposed to Solventy _foregsobutiont hminingys i Thus 
cross-connection of salt wells (and perhaps oil wells in 
future) by means of extending a hydrofracture from one well 
to the other has become a boon to the salt industry (Pullen 
Jr.,1958). The two wells so connected may be referred to as 
Source and target wells, with the pressurized well from 
which the hydrofracture originates as a source. The target 
well acts in one of two ways. Ideally the target well should 
act like a sink so that the propagating fracture will want 
to intersect it. In some cases however, with the high stress 
eoncentrationserarcundse@ity Che emaysactwasha barriermbotthe 
approaching fracture and most target wells act like minor 
sources, and thereby deflect the propagating fracture around 
themselves. To overcome this secondary role of the target 
well several techniques proposed have had promising effect. 
Htoewidlialways#abe jgqo00ds practice adtestoeroducecchydrauleae 
fractures from the target well either before or during the 
major fracture attempt from the source,(Bays et al, 1960). 
Back fracturing has proved to be successful in several cases 
where a fracture connection was not successful in the 
eriginalildmwvection. pit wsitalselgthought thartrmuchdofwthe 
stress concentration around the target well might be 
reduced, if the target well were pressured up and held 
during the fracture operation of the source well, at a 
pressure less than the breakdown pressure. Crawford and 


Collins (1954) point out the adverse effects on the 
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efficiency of a fluid injection project when injection and 
production wells are not oriented with the natural and 
induced fractures. If wells are properly oriented and 
fractures of adequate length are induced, radial flow is 
reduced and areal-sweep problems are minimised. 

We now consider the attempt to look critically at the 
lithologies in this study. The phenomenon of stress 
redistribution at depth in layered sequences of unlike 
rocks, known as arching, has been mentioned earlier. Arching 
may also be visualised in terms of massive competent beds 
acting as beams to carry the load imposed by crustal stress 
field and thereby reduce the stress on the less competent 
beds t aSuch@anierfectPisheimportant. to/Sthe Tdesigqnisots any 
Subsurface opening. Perhaps the greatest importance of 
arching is its influence through the stress distribution on 
the breakdown pressure. The object of a multiple fracturing 
process (multifrac) is to obtain greater fracture area 
through closer fracture spacing, than could be obtained from 
the same number of wells with single vertical fractures. 
This method requires Substantial deviation of the wells from 
the vertical and, therefore both the measured well length 
and drilling cost per unit length are increased when 
compared with a vertical well. Theory predicts a _ vertical 
hydraulically induced fracture in a strike-slip stress 
field, and such observations as are available confirm that 
in the Alberta sedimentary basin hydraulic fractures are 


generally vertical (Macleod, 1977 ). Thus it should be 
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possible to generate several fractures from a suitably 
deviated hole. Hydraulic fracturing theory predicts that 
fractures will propagate in a direction normal to the least 
principal ings Et stress or parallel to the maximum 
principal stress. Some experiments with non penetrating 
fluid indicate fractures may at least initiate parallel to 
the axis of the wellbore. A fracture initiated from an 
inclined open-hole wellbore mMighth ithene bresult ssintsa 
configuration unsuitable for the development of multiple 
vertical fractures. However, by casing the wellbore and 
perforating the caSing over a small depth range, .fracture 
initiagione Ganwebes imadearttorececur sfnom atpornt cortacsmahl 
Sphe rica hucavi bym(Steubhar cet baly.Al975). The orientation 
Should then be governed by the in-situ stress condition. 
Strubhar et al have shown experimentally that if the 
fractures are indeed oriented as dictated by in-situ 
stresses, the fractures will be sufficiently parallel and 
therefore be effective. The productivity index they found 
was two and three times that of a conventionally fractured 
well in the field. 
Conclusions 

1 has been argued earlier in this chapter and 
elsewhere that the consistency of the breakout azimuths 
irrespective of depth, lithology and age of the formation in 
almost vabliothe gwellsad and sethecriacta(tthataetthee breakout 
azimuths reflect only one out of the two joint systems of 


Babcock, lead to a safe conclusion that the breakouts are 
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stress controlled or possibly both. 
It 1S probable that some of the deviations from the 
generally northwest-southeast orientation of the breakouts 
can be attributed to local stress anomalies and some also to 
acerdentaliiydraubvcetensute Lnacturingeduringedgil linggwith 
high-density mud. Zoback and Pollard (1978) noted that 
although high viscosity fluid reduces fluid penetration into 
pre-existing fractureS or into permeable rocks, it has the 
disadvantage that such fluid may lead to hydraulic fractures 
Significantly in advance of the observed breakdown in the 
pressure-time history. 
On the hypothesis that breakouts are stress controlled, 
the following conclusions can be made. 
iIg@athetstaregonsstresstinacthe “earthisancrustrp iss. not ein 
general lithostatic. There are high horizontal stresses 
relative to the vertical stress irrespective of depth. 
Tectonic conditions and the various sedimentary layers 
may, however, influence this state of stress. 

2. Perhaps, some horizontal fractures were wrongly taken 
for vertical ones. Such horizontal fractures are 


believed to result from local inhomogeneities’ like 


bedding planes, joint-sets and the phenomenon of 
arching. 
3. Although breakouts tend to occur in preferential 


directions dictated by the orientation of the regional 
Stresses, induced vertical fractures may also be 


expected to have a preferential direction. In most cases 
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this direction 1s B@noesvdenctical to that of breakouts, 

but at right angles to them. 
Concluding remarks 

Therconsistency’ Of the datayin «thiss.region of North 
America and elsewhere on the continent strongly suggests 
that this meaSurement does not merely reflect local stresses 
but does, in fact, reflect a regional pattern. The existence 
of high stresses is not confined to any member of the rock 
formations, nor to any particular age of the formations. The 
phenomena of high stresses exist elsewhere in the North 
American continent and the world (Sbar and Sykes, 1973; 
Hastyet9/S 7174s sor) (OW 6 saMeGadretanduGay ,19¢7eseGough!,.21978- 
Zoback and Zoback, 1980) as evidenced by features like 
breakouts, folds, postglacial faults,and pop-ups. These 
stresses are probably related to the current tractions on 
the North American plate in the overall context of plate 


rectonics. 
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APPENDIX A 


A. Properties Of @ 
LS ial, ea O,' be the angles obtained from 


Zi,.eee+-,%, when the new zero direction is a. Let 


Cc’ = 1/n . Ecosg." ; S* = 1/n . Esin@, '------------------ A-1 
We have @.' = (@.-a)mod 27. So that 
C' = R.cos(G-a) ; S' = R.sin(@-a)--------------------- A-2. 


If we write (A-2) as 
C' = R'.cos@' ;: § = R'.sin@'--------------------------- A-3 
then we have 

@' = (@-a) mod 27 ; R'! = R --------------------------- A-4., 
Hence Equation 3-2 is satisfied. Also from (3-4) and (3-5) 
we have 

BN ee cr ee ea A-5 
which corresponds to the equation in the linear case 3(X. -X) 
= 0, 1.e.the sum of deviation about the mean vanishes. 
Circular Standard Deviation 

We have seen that S$, lies in the range (0, 1). To 

relate S, to the standard deviation oon the line, it is 
natural to use the following result for the wrapped normal 


distribution. 


iS jae oF, We can then define the circular standard 
deviation as 

a= 4 ee oe ie © Ree ere eee ei A-6 
The range of Vs of course (0,=). For «small S,, Equation 


A-6 reduces to g= (2.8,)”. a transformation first suggested 
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